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-Pace 

In  my  search  for  a  thesis  topic,  I  wanted  one  that  paralleled  my 
graduate  courses  in  the  flight  control  sequence.  The  topic  of  recon- 
figurable  control  laws  was  a  perfect  area  for  research.  This  is  a  new 
area  of  study  and  meets  the  needs  of  future  Air  Force  aircraft.  This 
study  was  proposed  by  the  AFWAL/rlight  Dynamics  Laboratory  (FDL)  and 
investigates  how  the  inherently  redundant  control  surfaces  of  future 
aircraft  could  be  utilised  after  a  surface  failure  by  employing  a 
digital  flight  control  system  (DFCS). 

This  thesis  provides  an  accurate  aircraft  model  for  further  re¬ 
search  into  the  area  of  reconfigurable  control  laws.  I  also  attempt  to 
show  how  a  direct  digital  design  method  using  the  entire  eigenstructure 
assignment  can  be  implemented  to  design  a  state  feedback  control  law. 

I  want  to  express  my  gratitude  to  Dr.  J.  D'Azzo  for  his  guidance 
in  the  thesis  development  and  unending  effort  in  the  reading  of  this 
thesis  for  accuracy  and  completeness.  My  appreciation  to  Mr.  Jerry 
Jenkins,  Flight  Dynamics  Laboratory  (FDL),  for  lis  assistance  in 
deriving  the  essential  aircraft  data  that  was  needed  to  complete  this 
report . 

I  also  wish  to  thank  my  sponsor,  Mr .  Duma  Ruboi  tus  (FDL), 
Professor  C.  H.  Houpis,  and  Capt.  J.  Sil\erthorn  for  their  assi'tunce 
and  suggestions. 
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This  thesis  investigates  control  of  an  aircraft  when  there  is  a 
primary  control  surface  failure.  The  object  of  this  study  is  to  recon¬ 
figure  the  remaining  control  surfaces  to  compensate  for  the  additional 
forces  and  moments  generated  by  the  inoperative  control  surface.  To 
study  this  flight  control  problem,  a  comprehensive  aircraft  model  is 
required  which  considers  each  control  surface  operating  individually. 

A  six  degree-of-freedom  aircraft  model  is  developed,  including  all 
the  individual  control  surfaces.  A  control  surface  input  can  produce 
both  a  lateral  and/or  a  longitudinal  response.  Thus,  the  equations  of 
motion  cannot  be  decoupled  for  the  design  of  the  control  laws.  The 
coupling  between  the  axes  requires  the  derivation  of  several  new  non- 
dimensional  control  derivatives.  Using  the  geometrical  properties  of 
the  aircraft  and  the  Digital  Datcom  computer  program,  the  needed  con¬ 
trol  derivatives  are  derived... 

With  a  comprehensive  aircraft  model  now  available-^he  entire 
eigenstructure  assignment  method  is  used  to  assign  both  the  eigenvalues 
and  the  eigenvectors  to  the  closed-loop  plant  matrix.  This  method  is 
used  for  the  dii'ect  digital  design  of  a  multivariable  discrete  regula¬ 


tor  and  tracker  control  law.  The  effect  of  increasing  the  number  of 
control  inputs  on  the  relative  degree  of  controllability  of  the  states 
was  determined  by  singular  value  decomposition.^ 

This  thesis  concludes  that  u  direct  digital  design  for  reconfigur¬ 
ing  the  multivariable  control  law  is  feasible.  However,  more  wind 
tunnel  data  is  essential  10  derive  the  additional  control  derivatives 
for  a  more  accurate  aix'craft  model  driven  by  individual  control 
surfaces.  Further  v/ork  is  also  necessary  to  perfect  the  assignment 
of  the  closed-loop  eigenvalues  and  eigenvectors. 
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DIRECT  DIGITAL  DESIGN  METHOD  FOR  RKC0NFIGURA3LE  MULTIVARIABLE 
CONTROL  LAWS  FOR  THE  A-7D  DIGITAC  II  AIRCRAFT 

I.  Introduction 

The  increased  flexibility  and  miniaturisation  of  digital  computers 
make  it  possible  to  use  digital  flight  control  systems  (DFCS)  on  future 
aircraft.  These  aircraft,  whether  they  are  inherently  stable  or 
unstable,  will  have  the  normal  sets  of  primary  control  surfaces 
(ailerons,  spoilers,  flaps,  and  horizontal  stabilizers),  split  into 
independently  controllable  surfaces.  The  combination  of  a  DFCS  and 
independently  movable  control  surfaces  suggest  that  a  method  for  rede¬ 
signing  the  control  laws  be  developed.  A  set  of  reconfigured  multi- 
variable  control  laws  could  then  be  stored  in  the  DFCS.  In  the  event  a 
primary  control  surface  becomes  inoperative  or  blown  away  due  to  battle 
damage,  a  reconfigured  multivariable  control  law  would  be  implemented 
by  the  DFCS  upon  the  detection  and  isolation  of  the  failure.  The 
design  method  developed  in  this  thesis  results  in  reconfigured  multi- 
variable  control  laws  using  state  variable  feedback  and  complete 
eigenstructure  assignment. 

This  chapter  presents  the  background,  problem,  scope,  assumptions 
and  approach  associated  with  this  thesis. 

Background 

Due  to  recent  advancements  in  digital  technology,  there  is  improved 
computational  speed  and  lower  cost  for  digital  computers  used  in  flight 
control  systems.  By  replacing  the  analog  flight  control  system  (AFCS) 
with  a  DFCS,  it  has  been  shown  one  can  achieve  a  faster  reaction  to  an 
input  disturbance  and  increased  flexibility  in  control  surface 
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utilisation  (Ref  10). 

The  DFCS  used  in  the  A-7D  Digitac  II  aircraft  is  an  experimental 
test  system  which  is  being  used  to  determine  if  the  mission  capabili¬ 
ties  are  improved  when  compared  to  the  AFCS.  The  test  program  is  also 
intended  to  provide  a  technology  base  for  future  development  of  digital 
and  multimode  flight  control  functions.  The  results  of  a  92  hour 
flight  test  program  conducted  between  February  1975  to  March  1976 
showed  that  the  DFCS  had  the  same  desired  high  degree  of  reliability 
and  flexibility  as  the  AFCS  and  provided  a  significant  improvement  in 
the  aircraft  handling  qualities  (Ref  10:9).  Therefore,  a  DFCS  with  the 
proper  detection  and  isolation  of  tie  loss  or  failure  cf  the  primary 
control  surface  is  capable  of  responding  faster  to  the  implementation 
of  a  reconfigured  control  law. 

Reconfiguration  of  the  control  laws  implies  reconfiguration  of  the 
control  surface  utilization  to  compensate  for  a  surface  failure  by 
battle  damage  and  using  the  remaining  surfaces  to  control  the  aircraft. 
In  the  past,  this  has  been  done  manually  by  pilots  with  some  success, 
particularly  in  aircraft  with  a  high  degree  of  aerodynamic  stability. 

In  the  future,  however,  the  loss  or  partial  loss  of  a  control  surface 
could  result  in  an  unrecoverable  flight  condition  in  inherently 
unstable  aircraft.  The  pilot  would  not  have  tine  to  identify, 
analyze,  and  take  corrective  action  after  a  control  surface  failure. 
Even  in  fiy-by-v/ire  (FB'.V)  controlled  configured  vehicles  (CCV)  of  the 
future,  manual  control  without  automated  assistance  to  achieve 
effective  reconfiguration  does  not  appear  feasible  (Ref  3:114).  Once 
a  reconfigured  control  law  is  implemented  and  control  of  the  aircraft 
is  regained,  tine  capabilities  and  response  of  the  aircraft  may  be 


2 


degraded.  Even  if  the  degradation  of  the  flying  qualities  only  allows 
the  pilot  a  get  home  and  land  capability,  the  reconfiguration  of  the 
control  laws  will  have  served  its  purpose  of  saving  the  pilot  and 
aircraft. 

Since  modern  aircraft  designs  tend  to  have  more  movable  control 
surfaces,  this  suggests  that  more  inherent  redundancy  is  available  if 
these  control  surfaces  can  be  properly  used.  With  a  DFCS  and  utiliza¬ 
tion  of  the  inherent  redundancy  of  the  control  surfaces,  it  is  feasible 
to  design  reconfigurable  control  lav/s. 

Problem 

The  Flight  Dynamics  Laboratory  is  presently  investigating  the 
feasibility  of  reconfiguring  the  aircraft  flight  control  laws  in  the 
event  of  a  primary  surface  failure.  For  example,  when  the  right 
aileron  or  the  left  horizontal  stabilizer  become  inoperative,  addi¬ 
tional  forces  and  moments  are  generated  that  are  not  normally  present. 
The  aileron  commands  would  produce  unwanted  longitudinal  forces  and 
moments,  while  stabilator  commands  would  produce  Lateral-directional 
inputs.  The  object  of  this  thesis  is  to  use  the  remaining  control 
surfaces  to  compensate  for  the  inoperative  control  surface.  To  study 
this  flight  problem,  a  comprehensive  aircraft  model  is  required  to 
consider  each  control  surface  individually.  At  present,  data  is  not 
available  for  consideration  of  each  individual  control  surface  input 
in  each  equation  of  notion. 

Scope 

This  investigation  is  limited  to  proposing  a  method  for  redesign¬ 
ing  the  control  laws  of  an  aircraft  after  one  of  the  primary  control 
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surfaces  becomes  inoperative.  This  study  determines  which  new  dimen¬ 
sional  control  derivatives  are  required  in  the  aircraft  equations  of 
motion  in  order  to  ensure  the  coupling  of  the  longitudinal  and  lateral 
motion  between  the  axes.  This  is  accomplished  by  the  use  of  additional 
independent  control  surface  inputs.  Next,  the  entire  eigenstructure 
assignment  is  used  to  specify  a  method  for  redesigning  the  control  laws. 
Only  continuous  time  simulation  is  used  in  evaluating  the  system 
responses  when  checking  a  system  configuration  or  control  law. 

Assumptions 

A  number  of  assumptions  are  made  to  simplify  the  complexity  of 
this  problem.  First,  the  aircraft  is  assumed  to  make  small  perturba¬ 
tions  about  a  trimmed  straight  and  level  flight  condition.  Secondly, 
the  aircraft  equations  of  motion  developed  are  used  in  the  stability 
axis.  Thirdly,  the  coefficients  of  the  control  derivatives  are  nor¬ 
mally  for  a  set  of  control  surfaces,  i.e.,  ailerons,  flaps,  horizontal 
stabilizers,  but  only  half  the  value  of  the  control  derivative 
coefficient,  with  the  proper  sign  is  used  when  the  control  surface 
set  is  split  into  two  independent  control  surfaces,  e.g.,  left  and 
right  aileron. 

Approach  and  Presentation 

The  design  approach  or  method  for  reconfiguring  the  control  laws 
in  this  thesis  looks  at  using  all  control  surfaces  independently  of 
one  another.  This  means  that  each  aileron,  each  horizontal  stabilizer, 
each  spoiler,  and  each  flap  is  controlled  independently.  This  allows 
for  maximum  utilization  of  all  the  surfaces  when  a  primary  surface 
fails.  To  make  use  of  ihese  additional  control  surfaces,  tiie  lateral 
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and  longitudinal  inputs  are  present  in  all  the  state  equations  of  the 
aircraft  model.  By  doing  so,  one  can  compensate  for  a  surface  failure 
by  using  state  variable  feedback.  To  accomplish  this,  additional 
control  derivatives  are  required  to  realize  the  effects  of  the  addi¬ 
tional  control  inputs  in  each  of  the  axes  as  the  aircraft  equations 
are  developed.  With  the  additional  control  inputs,  the  degree  of 
controllability  of  the  (A,B)  matrix  pair  is  examined  to  determine  if  a 
state(s)  requires  a  large  magnitude  of  control  input  in  order  to 
retain  control  of  the  aircraft  after  a  surface  failure. 

With  a  controllable  and,  perhaps,  a  reduced  order  plant,  the 
reconfigured  control  lav/  is  designed  using  the  complete  eigenstructure 
assignment  (CESA)  design  technique.  The  interactive  program  CESA  is 
used  to  generate  the  state  feedback  design.  This  is  an  iterative 
process  of  assigning  the  closed-loop  eigenvalues  and  selecting  the 
eigenvectors  from  the  controllable  subspace  for  the  assigned  eigen¬ 
values.  Then,  a  continuous  simulation  analysis  of  the  closed-loop 
system  responses  checks  the  reconfigured  control  law  implemented. 

This  thesis  is  composed  of  five  chapters.  Chapter  II  develops  the 
aircraft  model  and  references  Appendices  A  for  the  derivation  of  the 
aircraft  equations  of  motion.  The  method  used  to  obtain  the  additional 
control  derivatives  is  described  in  the  Appendix  B.  Because  each 
segment  of  the  control  surface  is  operated  independently,  there  is 
coupling  between  the  lateral  and  longitudinal  axes.-  This  coupling  is 
not  encountered  when  the  two  ailerons,  the  two  horizontal  stabilizers, 
and  the  two  spoilers  are  controlled  from  the  same  control  signal. 
Chapter  III  presents  the  theory  associated  with  entire  eigenstructure 
assignment  for  distinct  and  multiple  eigenvalue  assignment.  A 


regulator  design  example  is  used  to  illustrate  that  the  assigned 
closed-loop  eigenvalues  and  eigenvectors  can  be  obtained  from  a 
generated  control  lav/.  Chapter  IV  looks  at  the  modifications  required 
for  the  interactive  program  CESA.  Also  discussed  are  the  problems 
encountered  with  the  controllability  of  the  (A,B)  matrix  pair  and  the 
solutions  used  to  obtain  a  higher  degree  of  controllability.  Chapter  V 
presents  the  conclusions  and  recommendations. 


t 


6 


II.  The  Aircraft  Model 


Introduction 

An  accurate  model  is  desirable  for  any  control  design  methodj 
however,  limitations  and  assumptions  are  necessary  to  derive  a  model 
that  is  manageable  and  accurate.  This  thesis  uses  a  six  degree-of- 
freedom  (D-O-F)  aircraft  model.  The  model  includes  the  effect  of  in¬ 
dividual  control  surfaces,  many  of  which  are  traditionally  considered 
to  act  together  as  a  set.  For  example,  this  thesis  considers  the 
right  and  left  horizontal  stabilizer  as  independent  control  surfaces. 

A  linear  model  is  formed  by  linearizing  the  nonlinear  6  D-O-F  equations 
of  motion  about  a  nominal  operating  point  or  flight  condition. 

This  chapter  presents  a  general  description  of  the  A-7D  aircraft, 
a  discussion  on  the  linearized  aircraft  equations  of  motion  and  a  dis¬ 
cussion  on  the  new  control  derivatives  required  for  the  aircraft  model. 

General  Description 
. —  *  '  . — '  ■ 

The  A-7D  is  a  single  seat,  land  based,  light  attack  aircraft  with 
moderately  swept  wing  and  tail  surfaces.  The  fuselage  is  45.4  feet 
long  with  an  underslung  nose  inlet  duct,  and  is  powered  by  a  single 
TF41-A-1  engine  (Ref  8:3.1).  The  wing  area  is  375  square  feet  and  the 
mean  geometric  chord  (m.g.c.)  is  10.8  feet.  The  horizontal  stabilizer 
is  a  slab  tail  with  56.2  square  feet  of  area  and  a  m.g.c.  of  6.1  feet. 
The  distance  from  25  percent  of  the  wing  m.g.c.  to  25  percent  of  the 
tail's  m.g.c.  is  16.2  feet.  At  a  cruise  configuration,  the  aircraft 
weight  is  25,238  pounds,  and  is  flying  at  0.6  Mach  at  an  altitude  of 
15,000  feet.  Details  of  this  flight  condition  are  given  in 


System  Model 


The  aircraft,  model  in  this  thesis  uses  the  6  D-O-F  equations  of 
motion  for  an  accurate  model  without  approximations.  The  development 
of  the  aircraft  equations  is  in  Appendix  A.  In  the  derivation  of  the 
aircraft  model,  a  method  is  presented  to  resolve  the  problem  of  how  to 
achieve  control  of  the  lateral-directional  motion  with  a  longitudinal 
control  surface  and  the  control  of  the  longitudinal  motion  with  a 
lateral  control  surface.  This  is  a  primary  requirement  in  order  to 
accomplish  the  reconfiguration  of  the  control  law  after  the  loss  of  a 
primary  control  surface. 

To  achieve  this  control,  the  dimensional  control  derivatives  not 
normally  associated  with  either  the  latei'al-directional  or  longitudi¬ 
nal  axes  must  be  considered.  This  coupling  between  the  axes  is 
realized  through  the  non-dimensional  control  derivatives  associated 
with  the  force  and  moment  equations  control  surface  Inputs  for  the  non- 
traditional  surfaces.  For  example,  Equation  (l)  is  a  lateral- 
directional  input  coupled  into  a  longitudinal  control  derivative 
equation  by  Cn  . 


Sq  C 


D 


m 


(1) 


The  idea  cf  r.e  abo""  -jxample  applies  to  ail  xhe  new  dimensional 
control  derivatives  in  the  lateral-directional  and  longitudinal  con¬ 
trol  surf.ee  inputs  of  the  aircraft  equations  of  motion  as  shown  in 
Appendix  A. 
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By  including  the  individual  control  surfaces,  the  equations  cannot 
be  decoupled  to  separately  describe  the  longitudinal  and  lateral- 
directional  motion.  One  now  requires  knowledge  of  the  non-dimensional 
control  derivatives  for  these  control  surfaces  for  which  data  is  not 
usually  available.  The  next  section  provides  a  method  for  deriving 
this  data. 


Control  Derivatives 

To  successfully  carry  out  the  coupling  between  the  axes  requires 

the  derivation  of  several  non-dimensional  conti'ol  derivatives  for  the 

lateral-directional  and  longitudinal  control  inputs  of  the  equations  of 

motion.  The  control  derivatives  that  are  not  given  in  (Ref  11)  or 

(Ref  1)  are  derived  in  Appendix  B.  These  new  control  derivatives  were 

derived  by  digital  DATCOM,  (Ref  14),  and  conventional  aerodynamic 

1  techniques.  For  example;  the  rudder  was  modeled  as  an  asymmetrical 

wing  and  was  derived  from  digital  DATCOM  data.  This  is  explained 
6 

r 

in  more  detail  in  Appendix  B. 

Summary 

In  this  chapter,  a  physical  description  of  the  A-7D  is  presented. 
Then,  a  development  of  the  aircraft  equations  of  motion  is  presented  to 
introduce  a  method  that  achieves  coupling  of  the  lateral  and  longitud¬ 
inal  motion  into  the  opposite  axes  by  using  all  the  individual  control 
surfaces  and  new  non-dimensional  control  derivatives  for  the  respective 
axes.  With  a  comprehensive  system  model,  the  next  step  is  to  px*esent  a 
method  to  design  a  reconfigurable  multivariable  control  law  using  the 
,  complete  eigerstructure  assignment  design  technique. 
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III.  Theory  of  Entire  Eigenstructure  Assignment 
Introduction 

The  design  method  in  this  thesis  for  obtaining  the  reconfigured 
control  laws  uses  the  technique  of  entire  eigenstructure  assignment 
for  multi-input  multi-output  (MIMO)  control  systems.  The  assignment  of 
the  entire  eigenstructure,  consisting  of  both  eigenvalues  and  eigen¬ 
vectors,  is  illustrated  by  application  to  a  system  design  in  which 
distinct  eigenvalues  are  assigned.  The  use  of  multiple  eigenvalue 
assignment  is  particularly  important  for  sampled-data  systems  since  it 
can  produce  closed-loop  system  responses  having  finite  settling  times. 
Fully  exploiting  the  properties  of  a  multivariable  system  make  it  a 
powerful  design  tool  for  the  synthesis  of  sampled-data  controllers. 

This  enables  the  synthesis  of  a  sampled-data  control  system  which  con¬ 
sists  of  a  continuous-time  plant  and  a  digital  controller  that  produces 
a  control  input  signal  which  is  piecewise-constant  for  each  sampling 
period . 

This  chapter  presents  a  comprehensive  review  of  the  above  design 
technique  which  assigns  both  the  eigenvalue  spectrum  and  eigenvectors 
associated  for  the  closed-loop  plant  matrix.  First,  the  principles  of 
the  entire  eigenstructure  are  developed  for  the  case  of  distinct 
eigenvalues  and  then  a  regulator  is  synthesized  to  illustrate  these 
principles.  Then,  the  principles  for  multiple  eigenvalue  assignment 
are  presented.  This  is  followed  by  applying  the  theory  to  systems  •• 
augmented  with  integrators  for  the  design  of  a  tracker  controller  and 
a  discussion  on  continuous  simulation  (Ref  6). 
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Distinct  Eigenvalue  Assignment 


The  linear  multivariable  time-invariant  vector  matrix  equations 
governing  the  plant  and  its  output  in  the  continuous -time  domain  have 
the  forms 


x(t)  =  Ax(t)  +  B  u(t) 

(2) 

y(t)  =  Cx(t) 

(3) 

and  is  illustrated  in  Figure  1.  To  form  the  sampled-data  system,  a 
linear  time  invariant  discrete-time  MIM)  system  is  described  by  the 
matrix  difference  equation  (Ref  6:491) 

x(kT+T)  »  F(T)  x(kT)  +  G(T)  u(kT)  (4) 

yOcT)  =  C  x(kT)  (5) 

where  T  is  the  sampling  time  and  is  shown  in  Figure  2.  The  control  law 
for  applying  state-variable  feedback  for  the  discrete-time  system  is 

u( k'f )  -  Kx(kT).  (6) 

The  purpose  in  applying  state  feedback  is  to  assign  a  closed-loop 
self-con jugate  distinct  eigenvalue  spectrum 

o(F*GK)  =  {At,  A2 , - ,  An}  (7) 

and  the  associated  set  of  eigenvectors 

(V  XA2 . \  1  '  (8> 

>  " 

The  eigenvalues  and  eigenvectors  for  the  closed-loop  system  are  related 
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i 

i 

t 


n 


by  the  equation 


(F+GK)  Xx  =  XiXjL  (i  =  1,  2,  — ,  n) 


(9) 


which  my  be  put  in  the  form 


[F-X±l,  G] 


L  AiJ 


»■  0  (i  ~  1,  2,  — ,  n) 


(10) 


where 


\  =  (u) 


T  T  T 

and  where  [X  ,  ui  ]  is  a  vector  which  lies  in  the  kernel  or  null 

A ,  A  . 


space  of  the  matrix: 


S(Xi)  =  [F-XiI,G]  (i  =  1>  2,  --->  n)  (12) 


It  follows  from  Equation  (11)  that  the  required  state  feedback  matrix  K 
is  given  by 


K  = 


2 


(13) 


The  matrix  K  *8  real  and  simultaneously  assigns  both  the  desired  set  of 
eigenvalues  and  the  assigned  set  of  associated  eigenvectors. 

Equation  (il)  may  also  be  solved  directly  without  the  requirement  to 
obtain  the  inverse  matrix  shown  in  Equation  (13)  (F.ef  22).  The  ker 
3( )  (i  -  1,  2,  — ,  n)  imposes  constraints  on  the  eigenvector(s) 
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that  may  be  associated  with  the  assigned  eigenvalue  Aj(i  =  1,  2,  — ,n) 

by  identifying  a  subspace  within  which  the  eigenvector 

X.  (i  =  1,  2,  — ,  n)  must  be  located.  The  eigenvectors  must  be 
Ai 

linearly  independent  and  self-con jugate  so  that  the  inverse  matrix  X"^ 
in  Equation  (13)  exists  and  is  real.  Within  these  constraints,  there 
is  freedom  in  assigning  the  eigenstructure  of  the  closed-loop  plant 
matrix  F+GK.  This  is  illustrated  with  the  following  example  (Ref  6:13). 


Example  II-l. 

A  continuous  system  is  represented  by  the  following  state  and  out¬ 
put  equations: 


'0,1,0  0,0 

x(t)  =  0,0,1  x(t)+  1  ,  0  u(t) 

-6  ,-11,  -6  0  ,  1 


t)  =  [1  ,  1  ,  Ol 

.0  ,  0  ,  1 J 


The  design  of  the  regulator  in  this  example  is  carried  out  with  the 
interactive  computer  program  CESA  (Ref  8).  The  discrete-time  domain 
representation  is  obtained  with  a  sampling  period  of  T  =  1  second. 

The  eigenvectors  lie  in  the  null  space  of  the  matrix  S(A)  of  Equation 

(12). 


0.7474- A,  0.4530  ,  0.0735  ,  0.3261,  0.0421 
j(\)  -  0.4410  ,-0.0611-A,  0.0121  ,  0.4530,  0.0735 
0.0723  ,-0.5735  ,  -0.1334-A,  -1.0610,  0.0121 


The  open- loop  i impled-duta  eigenvalue  spectrum  for  the  discrete 
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plant  matrix  F  is 


o(F)  =  {0.50,  0.135,  0.368}  .  (16) 

The  feedback  matrix  K  is  to  be  computed  so  that  the  closed-loop  system 
is  asymptotically  stable  and  has  the  eigenvalue  spectrum 

o(F+GK)  =  {0.05,  0.1,  0.15}  .  (17) 


The  size  of  the  null  space  generated  is  dependent  upon  the  column 
dimension  of  the  input  matrix.  In  this  example,  the  null  space 
dimension  is  2  as  is  shown  for  the  specified  eigenvalues. 


ker  S(0.05)  =  span 


ker  S(0.10)  =  span 


ker  3(0.15)  =  span 


/ 

I~  — 

0.11152 

-0.33380 

1 

0.00005 

0 

. 

—  — i 

< 

0.19090 

-0.43770 

1 

0.00364 

0 

. 

0 . 27612~ 
-0.52460 

1 

-0.00229 

0 

. 

L*.  -J 

> 


/ 


> 


s 


-0.34692 

0.68090 

0 

-0.33305 

1 


> 


-0.32759 

0.57450 

0 

-0.27684 

1 


/ 


> 


-0.32393 

0.50770 

0 

-0.24100 

1 


> 


(18) 


(19) 


(20) 


From  each  null  space,  one  eigenvector  is  selected  ensuring  that 
the  requirement  for  linear  independence  is  met  and  also  satisfying 
Equation  (10).  Thus,  making  the  permissible  assignment  from 
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Equations  (18),  (19),  and  (20)  of 


XX,  XX2  **3 


“Xi  %  wx3 


-0.34692 

0.68090 


0.19090 

-0.43770 


-0.3330b  ’  0.00364 


0.27612 

-0.52460 


-0.00229 


The  required  feedback  matrix  is  computed  from  Equation  (13)  as 


K  =  -0.33305  0.00364  -0.00229  -0.34692  0.19090  0.27612 


0.68090  -0.43770  -0.52460 


The  resulting  closed-loop  sampled-data  plant  matrix  is 


a  -1.186,  -1.094,  -0.248?] 

3.127,  3.062,  O.7432] 


FyGK  = 


0.4992, 

0.7437, 

1.2240, 


0.2253, 

-0.3315, 

0.6238, 


0.0237 

-0.0459 

0.1393 


The  closed-loop  plant  matrix  has  the  required  eigenvalue  spectrum  of 
Equation  (16)  and  the  assigned  set  of  eigenvectors  of  Equation  (21),  as 
required.  Thus,  the  method  of  entire  eigenstructure  assignment  permits 
the  eigenvalues  to  be  altered  or  to  be  left  unchanged,  as  required. 

Even  when  an  eigenvalue  is  left  unchanged,  the  associated  eigenvector 
may  be  reassigned  within  the  appropriate  subspace. 

The  continuous-time  responses  of  the  states  for  the  closed-loop 


sampled-data  system  are  plotted  in  Figure  3  for  the  initial  conditions 


1 


x(o)  = 


1 

-1.5 

1 


(24) 


L 


The  plant  states  and  outputs  essentially  reach  steady-state  in  approxi¬ 
mately  three  sampling  periods. 


Multiple  Eigenvalue  Assignment 

The  assignment  of  multiple  eigenvalues  is  especially  useful  for 
sainpled-data  systems  since  the  assignment  of  all  eigenvalues  to  the 
origin,  i.e.,  X^  =  0(i  =  1,  2,  — n),  results  in  time-optimal  or  finite 
settling  time  responses.  When  it  is  desired  to  assign  multiple  eigen¬ 
values,  the  associated  eigenvectors  are  related  by  the  equation 


(F+G)Xi  -  xj,Xi  (i  =  1,  2,  —  p<n ) .  (25) 


Where  p  is  the  number  of  different  eigenvalues  and,  therefore,  only 
linearly  independent  eigenvectors  may  be  generated  from  the  null  space 
(Ref  6:2.23,  25).  Depending  upon  the  desired  algebraic  multiplicity, 
M^<n,  and  the  number  of  subblocks  desired  in  the  Jordan  normal  form  of 
the  closed-loop  matrix  for  each  assigned  eigenvalue  X^,  it  may  be 
necessary  to  generate  chains  of  generalized  eigenvectors  using 
the  equation 


(F+GK)  X. 


i+k 


=  X.x.  ,  +x.  , 

i  i+k  i+x- 


(k  --  1, 


— ,  Mi-1) 


(26) 


SECONDS ) 


(SECONDS) 


REGULATOR  DESIGN  EXAMPLE 


SECONDS) 


SECONOS) 


This  equation  may  be  put  in  the  form 


[F-XiI,G] 


■*» 

x  (mji» J ) 
Ai 

w.  jl 
Ai 


-1J) 


(J  ~  i >  2,  -*■-*>  k^;  i  =  1,  2,  — ,  p). 


(27) 


This  generates  strings  of  vectors  associated  with  the  eigenvalue  X^, 

where  is  the  1th  vector  in  the  jth  string  which  is  of  length 

since  the  ker  S(Xi)  occupies  a  subspace  of  size  equal  to  m,  the 
number  of  controls. 

Equation  (26)  allows  the  generation  of  eigenvector/generalized 
eigenvector  chains  from  the  original  null  space  given  by  Equation  (12) 
so  that  the  size  of  the  largest  control  index,  dj,  of  the  open-loop 
system  may  be  retained  for  the  closed-loop  system.  This  is  possible 
if  (l)  the  eigenvalue  spectrum  of  the  closed-loop  system  is  assigned 
to  the  origin,  (2)  the  combined  length  of  the  first  eigenvector/ 
generalized  eigenvector  chain  is  set  equal  to  d^  In  this  case,  the 
outputs  are  equal  to  the  inputs  after  a  finite  number  of  sampling 
times,  that  is,  the  system  is  deadbeat.  This  allows  a  deadbeat 
response  to  take  place  in  q  sampling  times,  where  q  is  the  order  of 
the  minimal  polynomial  of  the  closed-loop  system  and  is  equal  to  the 
largest  control  index  of  pair  (F,G)  (Ref  7:19). 

The  foregoing  theory  of  assigning  eigenvalues  and  selection  of 
eigenvectors  to  the  design  of  a  regulator  is  further  applied  to  the 
tracker  design  by  augmentation  of  the  system  with  integrators. 
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Tracker 

The  synthesis  of  a  tracking  system  which  is  required  to  track  a 
command  input,  uses  proportional  plus  integral  state  feedback  in  multi- 
variable  control  systems.  This  assures  zero  steady-state  error  for 
step  inputs  for  the  tracking  system  if  the  augmented  system  remains 
controllable  (Ref  4).  The  controllable  multivariable  plant  governed 
by  discrete-time  matrix  difference  and  output  equations  are 

xCkT+T)  =  Fx(kT)  +  Gu(kT)  (28) 


and 


y(kT)  =  Cx(kT)  (29) 

The  system  is  augmented  by  the  introduction  of  a  discrete-time  vector 
integrator  and  comparator 

z( kT+T)  =  z(kT)  +  v(kT)  -y(kT)  (30) 

Where  the  pxl  output  vector  y(kT)  is  required  to  track  the  pxl  piece- 
wise  constant  command  input  vector  v(kT).  The  number  of  command  inputs 
is  less  than  or  equal  to  the  number  of  outputs.  Inserting  Equation  (29) 
into  (30)  yields 

z(kT+T)  -  -Cx(kT)  +  z(kT)  +  v(kT)  (31) 

Thus,  from  Equation  (28)  and  (31),  the  discrete  composite  open- 
loop  augmented  tracking  system  is 

u(kT )  +  I  °  v(kT) 

**  V  ( 32 ) 


*x(  kT+T)’ 

*F  ,  0  ' 

x(kT) 

'g' 

z{  kT+T ). 

.-c  ,  I 

z(kT) 

.0. 
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Where  I  is  a  pxp  identity  matrix  and  p  is  the  number  of  integrators 
P 

used. 

The  proportional  plus  integral  control  law  is 

u(kT)  =  Kx(kT)  +  Kjs(kT)  (33) 

where  K  is  a  mxn  matrix,  Kj  is  a  mxp  matrix,  and  m<n.  The  closed-loop 
tracking  system  is  represented  by 


* x( kT+T ) 

Pf+GK,  GKn  rx(kT)l  ( 

3 

v(kT) 

t 

. z( kT+T ) . 

L-C  ,  IJ  [z(kT)J 

(34) 

and 

y(kT)  =[0,0] 

The  output  matrix  (C  ,  0]  in  Equation  (35)  is  augmented  with  zeros  for 
z(kT),  since  integrator  outputs  (Ref  6)  are  not  system  outputs. 

When  K  and  are  generated  that  produce  an  asymptotically  stable 
system,  the  closed-loop  eigenvalues  are  within  the  unit  circle.  The 
integrator  state,  then  reach  constant  values  with  a  constant  command 
v(kl'),  therefore, 


x(kT) 
I  z(kT) 


(35) 


lim  (z(kT+T )  -z(kT)]  =0  (36) 

k-x» 


and  from  Equation  (30), 

lim  [y( kT)  -v(kT)]  =  0 

k-x»  (37) 
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Therefore,  the  desired  state  outputs  will  track  the  command  inputs 
with  aero  steady-state  error. 

Continuous  Simulation 

Since  the  aircraft  plant  is  a  continuous  system,  it  is  necessary 
to  perform  a  continuous  simulation  of  the  closed-loop  system  in  order 
to  evaluate  the  effectiveness  of  the  digital  control  law.  This  is 
necessary  because  the  complete  closed-loop  system  is  partly  discrete 
and  partly  continuous.  Therefore,  the  system  is  simulated  using  the 
continuous  state  space  representation  with  piecewise-constant  control 
inputs . 

In  order  to  determine  the  continuous-time  solution  to  the  state 
space  iepresentation  of  the  linear  time  invariant  system,  a  discrete 
approximation  to  the  exact  solution  is  used.  This  method  is  used  to 
discretize  the  plant,  as  explained  in  the  preceding  section  on  sampled 
data  transformation,  with  a  sampling  time  that  is  less  than  or  equal 
to  TSAMP/12,  where  TSAMP  is  the  sampling  time  of  the  digital  computer 
used  in  the  control  law  (Ref  7:24).  The  state  and  output  responses 
are  obtained  from  one  sampling  time  to  the  next  with  the  control  input 
u(kT)  held  constant.  Then,  the  control  input  is  updated  at  the  next 
sampling  time  and  the  process  is  continued.  Therefore,  it  is  possible 
to  determine  t.he  system  responses  between  the  sampling  times.  This 
simulation  is  incorporated  in  the  computer  program  CESA  (Ref  8). 

Summary 

The  use  of  the  entire  eigeos tructure  assignment  design  technique 
for  multivariable  control  systems  is  possible  by  specifying  the 
desired  entire  eigenstructure  of  a  controllable  closed-loop  system. 
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There  is  complete  freedom  in  assigning  the  eigenvalues  and  the 
corresponding  eigenvectors  must  be  assigned  v/ithin  a  specified  null 
space.  Then,  a  state-variable  feedback  control  law  is  determined  for 
the  discrete-time  domain  system. 

This  design  method  can  be  applied  to  regulators  where  it  is 
desired  to  bring  the  system  back  to  the  steady  state  condition.  Also, 
the  system  oan  be  augmented  with  integrators  for  a  tracker  design 
where  the  output  tracks  a  command  input  vector  with  zero  steady-state 


error. 


IV.  Multivariable  Control  Law  Design 


Introduction 

The  reconfiguration  of  the  control  laws  for  an  aircraft  must  com¬ 
pensate  for  the  forces  and  moments  generated  by  a  partial  primary 
control  surface  failure.  The  design  was  achieved  by  using  the  complete 
eigenstructure  assignment  design  technique  for  a  multivariable  system. 
First,  it  was  necessary  to  derive  the  new  control  derivatives  for  the 
additional  independently  variable  control  surface  inputs  in  the  six 
degree-of-frecdom  (DOF)  aircraft  equations  of  motion.  In  accomplishing 
the  design,  not  only  must  the  controllability  of  the  (A,B)  matrix  pair 
be  determined,  but  the  degree  of  controllability  is  also  important. 

The  degree  of  controllability  for  this  multivariable  system  is  impor¬ 
tant  since  this  determines  the  magnitude  of  the  feedback  matrix 
coefficients  and  the  magnitude  of  the  plecewise-constant  control  input 
signals.  Should  there  be  a  state  with  a  low  degree  of  controllability, 
it  is  possible  to  remove  that  state  from  the  model  representation  if 
the  response  o°  the  original  system  is  not  affected. 

This  chapter  presents  the  modification  that  was  required  to  make 
the  interactive  program  CESA  compatible  with  the  control  input  matrix 
o1'  the  aircraft  model.  Then  the  degree  of  controllability  of  the 
system  Is  examined  by  using  a  method  that  balances  the  system  matrices 
in  a  new  state  coordinate  system,  allowing  for  state  reduction,  if 
desired.  Having  established  the  degree  of  controllability  of  the  • 
model  which  was  developed,  the  design  scheme  for  a  multivariable 
discrete- time  tracker  and  regulator  is  presented. 


CESA  Modification 


The  interactive  computer  program  CESA  provides  the  user  with  an 
interactive  tool  to  design  the  control  laws  for  regulators,  disturb¬ 
ance  rejectors,  and  trackers.  The  program  has  the  capability  for 
designing  a  multivariable  control  law  and  then  evaluating  the  system 
responses  by  continuous  or  discrete  simulation.  Should  the  control 
law  be  unacceptable,  it  can  be  redesigned  without  having  to  reenter 
the  system  matrices  and  sampling  times.  However,  the  closed-loop 
eigenvalue  spectrum  and  eigenvector  assignment  must  be  respecified. 
Thus,  the  system  data  is  stored  in  memory  for  an  intei*active  design 
process  (Ref  8:157). 

A  program  modification  was  required  due  to  the  dimension  limita¬ 
tions  on  the  continuous  and  discrete-time  control  input  matrices. 

Also,  other  matrices  and  vectors  used  in  the  design  process  in  the 

program  are  dependent  upon  the  number  of  control  inputs.  This  thesis 

required  that  the  control  inp  it  matrix  dimension  be  increased  in  size 
to  an  nxm  matrix,  where  ra  is  eight  and  this  does  not  include  the  failed 
surface.  This  allows  for  the  use  of  the  left  and  right  control  sur¬ 
face  input  of  the  unit  horizontal  tail  (UKT),  ailerons,  spoilers, 
flaps  and  the  single  vertical  tail  rudder.  To  accomplish  the  modifi¬ 
cation,  a  thorough  understanding  of  the  theory,  computer  program 
overlay  structure,  the  programs  and  subroutines  used  in  the  design, 
and  the  simulation  of  the  discrete  state  feedback  control  law 
generated  for  the  A-7I)  aircraft  model  was  necessary. 

Once  the  arguments  of  the  arrays  and  vectors  were  updated,  the 

program  was;  recompiled  for  use  on  the  AOD  Cyber  computer'  at  the'  Air 

Force  .’/light  Aeronautical  Laboratories  (AF.VAL)  using  the  following 
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external  libraries:  IMSL,  EISPACK,  CCAUX,  and  CCPLOT56X.  To  verify 

i 

and  validate  the  updated  version,  the  third-order  system  example  of 
Chapter  III  and  the  aircraft  model  of  Appendix  A  were  used  during 
testing.  Then,  if  further  corrections  or  parameter  changes  were 
necessary,  it  became  an  iterative  process  until  all  necessary  changes 
were  made  and  validated  to  complete  the  program  modification. 

Once  the  interactive  program  modifications  were  completed,  the 
effect  of  increasing  the  number  of  control  inputs  on  the  degree  of 
controllability  of  the  (A,B)  matrix  pair  was  investigated. 

Controllability 

It  is  reemphasized  that  the  interactive  program  CESA  checks  the 
controllability  of  the  (F,G)  matrix  pair,  but  not  the  degree  of 
controllability  of  the  plant  states.  One  is  interested  in  the  degree 
of  controllability  when  using  state  feedback  because  it  does  affect 
the  overall  gain  of  the  state  feedback  matrix.  If  the  state  feedback 
matrix  produces  a  multivariable  control  lav/  that  applies  large  control 
inputs  to  the  servo  actuators,  then  the  physical  limits  of  the  actua¬ 
tors  are  exceeded.  Therefore,  removing  a  state  of  low  controllability 
from  the  aircraft  model  improves  ( reduces )  the  average  gain  of  the 
feedback  matrix.  To  improve  the  controllability  of  the  model,  two 
approaches  are  available.  One  is  to  determine  the  magnitude  of  the 
singular  values  of  the  matrix  that  are  generated  by  the  Lfoore 
balancing  algorithm  (Ref  12:18-20).  This  has  been  implemented  in  the 
AFIT  program  MIMO  (Ref  13).  The  second  approach  is  to  reduce  the 
number  of1  control  inputs  to  only  those  control  surface  inputs  that 
have  the  greatest  effect  on  the  aircraft  response. 
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The  Moore  algorithm  (Ref  15,  16)  defines  a  state  space  coordinate 
system  which  orders  the  state  variables  with  respect  to  the  controlla¬ 
bility  and  observability  properties,  i.e.,  the  most  to  the  least 
controllable,  observable  state.  This  is  done  by  transforming  the 
system  into  an  internally  balanced  state  coordinate  system,  where 

x  =  Tx'  (38) 

Using  the  linear  transformation  matrix  (Ref  16) 


T 


-1  b 

u  i  u..£ 

oo  H  H 


(39) 


where  U 


H 


INF 


=  left  singular  vector  of  the  observability  grammian 
=  diagonal  matrix  •ontnining  singular  values  of  the 
observability  grammian 
*  left  singular  vectors  of  the  11^,,  matrix 

=  Eujui 

0  o  c  c 

-  diagonal  matrix  containing  singular  values  of  the 

matrix 

INF 


Uq  =  left  singular  vectors  of  the  controllability  grammian 
1 1  diagonal  matrix  containing  singular  values  of  the 

G 

controllability  grammian 

This  transformation  to  the  internally  balanced  state  coordinate  system 
results  in  the  following  equations: 


x'  (t)  -  A'xtt)  ♦  B'u(t) 

y(t)  =  C'x'(t)  (40) 
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where  A"'=T  AT,  B“  =  T  B,  and  C'  =  CT  .  However,  to  retain  the 
physical  inputs  and  output  positions  of  the  original  system,  which  are 
the  important  physically  measurable  quantities  desired,  the  B'  input 
matrix  is  premultiplied  by  '  and  the  C'  output  matrix  is  post  multi¬ 
plied  by  T”\  so  that  the  B'  input  and  C'  output  matrices  are  ordered 
as  in  the  original  system.  Then,  the  state  and  output  equations  are, 

x'(t)  =  A'x1t)+  B'u(t) 

y(t)  =  x'(t)  (41) 

Since  the  internally  balanced  state  coordinate  system  orders  the 
states  with  respect  to  the  controllability  and  observability  proper¬ 
ties,  then  deletion  of  the  bottom  state(s)  of  the  balanced  matrix 
representation  is,  in  effect,  stripping  away  the  least  controllable/ 
observable  redundant  state(s).  The  following  is  an  example  of  state 
reduction  from  a  tiurd-order  to  an  approximately  equivalent  second- 
order  system: 


where  A',  B"  ,  and  C'  are  the  reduced  matrices.  This  does  not  imply 
11  °R 

that  model  reduction  is  necessary  to  use,  the  internally  balanced 
matrices  generated  by  Moore's  algorithm  effectively.  The  balanced 
matrices  can  be  used  us  generated,  the  advantage  being  that  the 
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mathematical  properties  of  the  plant  matrix  are  improved.  An  example 
is  that  the  condition  number  of  the  aircraft  model's  plant  matrix  is 
reduced  from  140653  to  127  when  in  the  balanced  form. 

A  method  has  been  presented  which  determines  the  degree  of  con¬ 
trollability  of  the  states  and  permits  model  reduction  of  the  original 
system,  if  desired.  The  next  step  in  the  design  method  for  reconfig- 
urable  multivariable  control  law3  is  to  develop  a  state  feedback  matrix 
for  the  aircraft  model  of  Appendix  A  for  a  tracker  and  regulator  con¬ 
troller  using  the  entire  eigenstructure  assignment  design  technique. 

Control  Law  Design 

A  linear  multivariable  discrete-time  tracker  or  regulator  con¬ 
troller  is  designed  to  reconfigure  the  control  surface  response  to 
compensate  for  the  failure  of  a  primary  control  surface.  This  section 
discusses  the  problems  and  solutions  encountered  with  the  method  for 
the  tracker  and  regulator  controller  design. 

The  tracker  design  is  desired  because  the  output  tracks  a 
commanded  input.  With  a  primary  control  surface  failure,  once  it  is 
detected  and  isolated  by  a  fault  detection  system,  a  control  law  is 
implemented  from  a  family  of  control  laws  in  the  digital  flight  con¬ 
trol  system  (DFGS)  to  reconfigure  the  response  of  the  remaining 
independently  operating  control  surfaces.  The  reconfigured  tracker 
control  law  would  then  follow  the  command  inputs,  removing  the  unde¬ 
sired  response  caused  by  the  failed  surface.  The  open-loop  Equation 
(32)  for  the  tracker  design  is  composed  of  Equations  (A-22,  A-23,  A-24, 
and  A-25)  augmented  with  the  integrator  and  comparator  of  Equation  (31). 
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Figure  4.  Control  Surface  Failure  Vs.  Reconfigured 
Input  (9)  Matrix 

Since  the  selected  output  was  required  to  track  four  command  inputs, 
four  integrators  were  used.  This  design  can  be  used  for  the  failure  of 
any  primary  control  surface  by  reconfiguring  the  control  law  to  use  a 
set  of  control  inputs  as  shown  in  Figure  4. 

The  desired  tracker  control  law  is 

u(kT)  -  Kx(kT)  +  z(kT)  (44) 

where  the  feedback  matrix  [K,  Kj]  is  computed  so  that  the  closed-loop 
system  is  asymptotically  stable  and  has  the  selected  eigenvalue 
spectrum  for  the  tracker  containing  4  integrators,  see  (A-25), 

a[F  +  G  (K,  Ki )]  =  (0.01,  0.03,  0.05,  0.06,  0.07,  0.08, 

0.1,  0.15,  0.2,  r  i,  0.4,  0.6}.  (45) 

The  state  feedback  matrix  is  related  to  the  eigenvalues  and  associated 
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eigenvectors  of  the  closed-loop  system  by 


i 


♦ 


L  ^  » 


H 


r. 


I?  3 


[F  -  AiI,G] 


<u. 


_  0  (i  =  1,  2,  ■  ■ >)  n) 


where 


0)^  =  K 


(i  =  1,  2,  ...,  n) 


and 


(i  =  1,  2,  . . .,  n) 


is  a  vector  that  lies  in  the  null  space  of  the  matrix 


S  (A.)  =  [F  -A. > G]  (i  =  1,  2,  ...,  n) 


Then  the  feedback  matrix  is  given  as  (Ref  6:17) 


K  -  [«i,  oi2,  •••>  (Xi>  x2>  •••>  Xn] 


-1 


K  =  fix 


-1 


(46) 


(47) 


(48) 


(49) 


(50) 


As  stated  in  Chapter  III,  the  eigenvectors  must  lie  within  the  ker 
S(  A^ )  of  the  specified  null  space.  Also,  the  columns  of  the  modal 
matrix  X  must  be  independent  for  the  matrix  X-^  of  Equation  (50)  to 
exist.  Also,  if  the  matrix  X  of  Equation  (50)  is  ill-conditioned,  an 
inverse  cannot  be  generated,  but  an  erroneous  inverse  will  be  produced, 
i.e.,  XX  ^  1  I.  This  is  a  computer  computational  problem. 

The  initial  tracker  design  on  the  interactive  program  CESA  did 
generate  a  control  law  that  was  formed  using  an  incorrect  X~\  That 
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tracker  control  law  had  a  very  high  average  gain  of  approximately  109 
and  the  closed-loop  specified  eigenvalues  were  not  contained  in  the 
characteristic  equation  of  the  closed-loop  matrix  [F+GK] .  This  was  a 
result  of  the  problem  areas  mentioned  previously:  (l)  degree  of 
controllability  of  the  (F,G)  matrix  pair,  and  (2)  an  ill-conditioned 
matrix  X  of  Equation  (50).  This  result  was  not  anticipated  and,  at 
first,  it  was  not  known  if  the  subroutines  and  programs  in  CESA,  the 
external  reference,  or  the  problems  stated  above  caused  the  high  gains 
and  non-reproduction  of  the  closed-loop  eigenvalues.  Extensive  testing 
was  done  on  the  CESA  subroutines  since  a  larger  number  of  control 
variables  were  being  used  which  called  for  an  even  larger  augmented 
tracking  matrix  of  Equation  (32).  Also,  the  external  reference 
libraries  were  tested  to  ensure  that  the  null  spaces,  eigenvalues,  and 
X-1  being  generated  were  correct.  During  the  testing  it  was  discovered 
that  other  available  library  routines  can  calculate  the  feedback 
matrix  K,  i.e.,  solving  for  the  K  matrix  for  the  linear  algebraic 
equation  having  the  form  XK  =  0.  The  CESA  program  and  the  external 
references  functioned  as  intended  and  it  was  determined  that  the 
difficulty  lay  in  the  ill-conditioned  X  matrix  of  Equation  (47). 

The  controllability  matrix  associated  with  the  (F,G)  matrix  pair 
must  have  full  rank.  A  necessary  condition  for  the  system  to  be 
controllable  is  that  the  rank  of  the  controllability  matrix  be  equal 
to  n,  the  dimension  of  the  plant  matrix,  where  the  controllability 
matrix  is  defined  as 

Mc  =  [G,  FG,  ...,  Fn_1G]  (51) 

However ,  the  degree  of  controllability  of  the  states  is  unknown.  To 
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drive  a  state  to  the  desired  steady  state  condition,  the  magnitude  of 
the  feedback  gain  is  dependent  upon  the  controllability  of  that  state. 
The  feedback  matrix  is  related  to  the  eigenvectors  which  determine  the 
magnitude  of  the  transient  responses.  The  controllability  problem  was 
further  investigated  by  considering  the  regulator  design.  The  regula- 
tor  does  not  use  integrators  to  augment  the  system  and  involves  working 
directly  with  the  plant  matrix  and  the  control  inputs  in  determining 
the  degree  of  controllability  of  the  (A,B)  matrix  pair. 

Regulator  Control  Law 

For  the  remainder  of  this  study,  the  right  aileron  is  the  surface 
failure,  r.ua  it  i3  assumed  to  be  locked  in  a  neutral  position.  Before 
the  regulator  design  was  implemented,  the  Moore  internally  balancing 
algorithm  was  applied  to  Equations  (A-22,  A-23,  and  A-26)  using  the 
interactive  program  MIMO  (Ref  13).  The  time  responses  of  the  original, 
balanced  and  reduced  balanced  system  were  compared  to  ensure  that  a 
true  representation  of  the  model  was  l’etained.  Once  this  was  estab¬ 
lished,  a  regulator  multi variable  control  law  was  designed  for  the 
three  systems.  Then  the  time-  responses  of  the  three  systems  were 
compared. 

The  degree  of  controllability  of  the  states  was  determined  from 
the  singular  values  of  the  B, Jf,  matrix  by  means  of  the  Moore  algorithm. 
A  guideline  to  model  reduction  is  the  magnitude  and  grouping  of  the 
.singular  values.  The  singul  a’  values  of  the  matrix  were  computed 
for  impulse  input  balancing  fter  a  finite  time  of  20  seconds.  A 
finite  time  was  used  due  to  the  spiral  divergence  pole  in  the  right 
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half  S-planc.  The  singular  values 


for  the  plant  matrix  are, 


0l  =  14.896 
o2  -  13.617 
o  3  =  10.569 
o4  s  9.300 
a 5  -  8.851 
a6  =  3.945 
a7  =  3.278 

a8  -  2.500  (52) 

These  singular  values  of  the  plant  indicate  that  at  least  one  state 
corresponding  to  oa  and  possibly  two  others,  o6  and  07  could  be  removed 
due  to  a  low  degree  of  controllability.  The  least  controllable  state 
is  represented  by  the  value  for  o8.  To  establish  if  the  balanced 
matrices  affect  the  time  response  of  the  system,  the  original  system 
was  compared  with  the  balanced  and  reduced  balanced  system.  As 
illustrated  in  Figure  5,  the  output  responses  of  each  system  are 
approximately  equal.  Recall  that  C'T~^  is  Cq  and  T  P*  is  Bq,  which 
insures  the  physical  outputs  and  inputs  remain  the  same  is  the 
original  system. 

The  eigenvalues  of  the  original  plant  are,  as  follows: 

Xj,2  r  -0.4148E-2  +_  j  Q.8209K-1  phugoid 
A 3, 4  --  -0.4358  +  j  2.040  dutch  roll 

A 5, 8  -  -0.7949  +_  j  2.848  short  period 
\y  -  0.4934E-1  spiral  divergence 

Ag  =  -2.878  roll  subsidence  (53) 

The  time  response  plots  for  the  original  system,  Figure  6,  indicate 
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Output  Signals  of  (a)  Original  System 
and  (o)  Internally  Balanced  System 
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that  some  form  of  compensation  or  feedback  control  is  necessary  to 
correct  for  the  unstable  spiral  divergence  mode  present  in  the  pitch 
and  roll  angle  responses.  The  pitch  response  of  Figure  6a  is  for  an 
impulse  of  0.1  rad  and  has  a  maximum  negative  pitch  angle  of  0.045  rad 
and  an  initial  settling  time  of  6  seconds  before  slowly  diverging  with 
a  constant  positive  pitch  angle.  The  roll  response  of  Figure  6c  has  an 
initial  settling  time  of  11  seconds  before  a  constant  roll  angle 
divergence  takes  over.  To  illustrate  that  the  horizontal  stabilizer 
UHT  does  function  on  the  model  as  left  and  right  independent  control 
surfaces,  Figures  6e  and  6f  show  that  the  left  stabilizer  produces  a 
positive  roll  and  the  right  stabilizer  a  negative  roll  with  a  0.1  rad 
impulse  input.  This  illustrates  the  coupling  between  the  axes  as  a 
result  of  longitudinal  inputs  into  the  lateral-directional  equation  of 
motion. 

After  balancing  the  system  state  equations  with  the  Moore  algo¬ 
rithm,  the  A",  B'q,  and  C"  matrices  are  as  shown  in  Equation  (A-25 
through  A-27)..  The  eigenvalues  for  the  A'  matrix  remain  the  same  as 
for  the  physical  plant  A  matrix  Equation  (53)  and  are  as  follows: 

Ai  2  '  -0.4148E-2  +  j  0.8209E-1 
X3  4  =  -0.4353  +  j  2.040 
*5,6  =  -0.7949  +  j  2.848 
X7  =  0.4934E-1 

X8  =  -2.878  (54) 

where  Xj  to  X8  correspond  to  the  same  characteristic  roots  of  the  ori¬ 
ginal  system.  Figure  7  illustrates  that  the  time  response  plots  of  the 
system  equations  with  internally  balanced  matrices  remained  the  same  as 
the  original  system.  The  yaw  a  igle  response  of  Figure  7f  for  a  0.07  rad 
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Figure  6b.  Pitch  Rate  For  0.1  Rad  Impulse 
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Figure  6d.  Roll  Rate  for  0.2  Rad  Impulse  6  Input 
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?igure  6f.  Roll  Angle  for  0.1  Rad  Impulse  6,  Input 


Figure  6g.  Sideslip  Angle  for  0.07  Rad  Impulse  5  Input 


impulse  input  results  in  the  some  form  of  response  as  for  a  0.2  rad 
impulse  input,  Figure  6h,  as  in  the  original  system  but  of  a  different 
magnitude  due  to  the  impulse  input.  The  difference  between  the  original 
and  balanced  matrices  is  that  the  states  of  the  system  are  ordered  by 
their  degree  of  controllability.  Thus,  the  redundant  states  can  be 
removed  as  long  as  the  system  characteristics  are  not  destroyed. 

The  A'  matrix  of  Equation  (A-25)  was  reduced  to  seven  states,  as 
illustrated  in  Equations  (42  and  43)  and  as  indicated  by  the  singular 
values  of  Equation  (52).  The  characteristic  eigenvalues  of  the  reduced 
internally  balanced  system  are: 

\x  2  3  -0.4142E-3  +  J  0.8225E-1  phugoid 
*3\  =  -0.4731  +  j  1.959  dutch  roll 

A 5  6  =  -0.8050  +.  j  2.825  short  period 

A 7  =  0.4935E-1  spiral  divergence  (55) 

The  eigenvalues  of  the  reduced  balanced  system  have  shifted  a  small 
amount  and  the  roll  subsidence  mode  is  no  longer  present.  The  reduced 
system  tine  response  plots  are  shown  in  Figure  8. 

The  pitch  and  pitch  rate  time  response  of  Figures  8a  and  8b  re¬ 
mained  the  same  as  for  the  original  system.  However,  tae  roll  angle  of 
Figure  So  shows  a  phase  change,  an  increase  in  magnitude  from  a  max  of 
0.024  radian  to  0.13  radian,  and  the  settling  time  decreased  by  0.5 
seconds.  These  changes  are  attributed  to  the  shifting  of  the  poles  and 

zeros  of  the  <f>/ 5  transfer  function.  In  Figure  8d,  the  roll  rate 
al 

settling  time  was  increased  from  approximately  6  sec  to  10  sec,  while 
the  magnitude  decreased  and  the  dutch  roll  oscillation  decreased.  The 
yaw  rate,  Figure  8f,  also  indicates  a  change  in  phase  and  magnitude. 


‘igure  7a.  Pitch  Angle  for  0.1  Rad  Impulse 


Figure  7b.  Pitch  Rate  for  0.1  Rad  Impulse 


Figure  7e.  Roll  Angle  for  0.2  Rad  Impulse  Input 
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Figure  7e.  Sideslip  Angle  for  0.07  Rad  Impulse  6  Input 


Figure  7f .  Yaw  Angle  for  0.07  Rad  Impulse 


Figure  8a.  Pitch  Angle  for  0.1  Rad  Impulse 


Figure  3c.  Roil  Angle  for  0.2  Rad  Impulse  5  Input 
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Figure  Sd.  Roll  Rate  for  0.2  Rad  Impulse  6  Input 


Figure  8e.  Sideslip  Angle  for  0.07  Rad  Impulse  S  Input 
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Figure  8f.  Yaw  Angle  for  0.07  Rad  Impulse  <5  Input 


Even  with  the  system  response  changes  pointed  out,  the  reduced  balanced 
system  still  retains  all  the  characteristics  of  the  aircraft  model. 
Therefore,  the  removal  of  a  state  with  a  low  degree  of  controllability 
does  not  affect  the  output  response  of  this  aircraft  model. 

Having  established  the  degree  of  controllability  of  the  states  by 
the  Moore  algorithm,  the  number  of  control  inputs  was  reduced  to  improve 
the  controllability  of  the  (A,  B)  matrix  pair.  A  reconfigured  multi- 
variable  control  law  for  a  regulator  controller  was  developed  for  the 
original,  balanced,  and  reduced  balanced  systems  using  a  sampling  time 
of  0.5  seconds. 

The  family  of  reconfigured  control  inputs  to  compensate  for  a 
failed  surface  in  the  regulator  design  is  shown  in  Figure  9. 
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Figure  9.  Control  Surface  Failure  Vs. 

Reconfigured  Input  (7)  Matrix 


A  block  diagram  for  the  state  variable  feedback  regulator  con- 
4»  troller  is  shown  in  Figure  10  in  which  all  the  variables  are  vector 

quantities. 

The  regulator  design  uses  Equations  (A-20)  and  (A-21)  except  for 
the  flap  inputs,  and  (A-24).  The  outputs  are  q,  0,  (5,  p,  r,  and  $. 

The  feedback  matrix  K  of  Equation  (47)  is  computed  so  that  the  closed- 
loop  original  system  is  asymptotically  stable  and  has  the  closed-loop 
eigenvalue  spectrum. 

a(F  ♦  GK)  =  {0.001,  0.002,  0.03,  0.3,  0.4,  0.625,  0.95,  0.99}  (56) 

The  vectors  assigned  from  the  null  space  of  Equation  (49)  for  each  X^ 
are: 


(57) 


•-.be  X  matrix  of  eigenvectors  is 


0. 

0. 

-0.1715E-03 

0. 

O.B531 

-0.4633 

0.4125 

0.7000 

0. 

1.000 

-0.5388 


1.000 

0. 

0. 

0  373E-01 

-0.31-2E-01 

0. 

0. 

0.1348E-01 

. 3257E-02 

-0.5443E-02 

-0.6069E-02 

392E-02 

0.2839 

-0.1677 

0. 

0 

0. 

1.000 

-0.3124 

6.1335E-03 

-0.2310 

0.1575 

0. 

0.3H7E-01 

0.1902E-01 

-0.1490E-01 

-0.7364E-03 

0. 

0.6591E-04 

0. 


0. 

-0.2891E-02 
-0.2430E-02 
.2315E-01 
-0.364 
1.000 
0. 3950E-02 
-Q.2633H-01 


0. 

0.1693E-02 

0.1535E-O2 

-0.7605E-01 

0.6055 

0. 

-0.2L27E-01 

1.000 


and  the  associated  il  matrix  is 


~ 0.4294 

0.4592 

-0.6622 

0.1458 

0.6545 

-0.3835 

0.5430 

-0.1698 

-0.1273 

0.9901E-01 

0.2618 

-0.9987 

-0.4618 

0.3591 

-0.5707 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1.000 

0.5741 

-0.3831E-01 

0.3030E-02 

-0.8611 

0. 

0. 

0. 

0.3593 

0.4353 

-Q.1491E-01 

-0.1055 

0.3299 

0.1503 

0.5032E-01 

0. 

0. 

0. 

1  .000 

0. 

0. 

0. 

0. 

0. 

0. 

• 
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The  size  of  the  null  space  and  the  row  dimension  of  the  Q  matrix  is 
dependent  upon  the  number  of  control  inputs.  The  computed  feedback 


malrix  K  of  Equation  (50)  for  the  original  system  was 


0.4714 

-0.7043 

-0.5807 

0.3532 

-0.3755 

-0.3578 

0.4705 

-0.1745 

0.1000 

-0.1523 

0.2659E-G1 

0.7001 

0.3585 

0.1414 

-1.237 

-3.476 

-0.162/ 

15.98 

-6.394 

-42.42 

-0.1064 

2.852 

4.820 

12.54 

-0.1556E-01 

-0. 29572-01 

-0.9141 

-0.8416" 

0.5861 

0.95422-01 

0.1055 

0.5018 

0.3425 

-0. 71422-01 

0.2144 

0.1805 

1.329 

0.2368 

-1.252 

-1.095 

-1.199 

0.7815 

-0.9945E-01 

-2.519 

-8.155 

-1.241 

4.150 

5.967 

with  an  average  gain  of  2.97.  The  resulting  closed-loop  sampled-data 
plant  n.atrix  is 


~0.2000E-02 

-0.1705E-12 

0.2842E-13 

-0.9337E-02 

0.5041 

0.1754 

0.1018E-02 

0.9657E-01 

-0.2298 

-0.1351E-02 

-0.6097E-01 

0.2059 

0.6413E-03 

0.1243E-01 

0.1492E-01 

0.3197E-13 

0.1776E-13 

0.7105E-13 

0.2198E-02 

-0. 3059E-01 

-0.6962E-01 

0.5137E-03 

-0.9197E-02 

-0.1965E-01 

0.4547E-12 
-0.2796E-01 
-0.8865 
0.7611 
0. 3888E-01 
0.3553E-12 
-0.1904 
-0.5355E-01 


0.2842E-13 
-0.1748E-01 
-0.4761E-01 
-0.1562E-01 
0.7357 
-0.1990E-12 
0.6788 
0.2030  . 


-0.1065E-13 
-0.2436E-02 
0.1046E-01 
0.2336E-02 
-0.2537E-01 
0.9500 
0, 8965E-0! 
0.2595E-01 


-0.2842E-13 

0.2227E-01 

0.1352E-02 

0.1241E-02 

-0.1566 


-0.2919 

0.1769 


-0.9948E-13 

0.9482E-02 

-0.3685E-01 

-0.8076E-02 

0.1536 


0.1279E-12  0.1181E-12 


-0.4526 

0.8668 


The  closed-loop  plant  matrix  has  the  specified  eigenvalue  spectrum  of 
Equation  ( 58 ) . 

A  continuous-time  simulation  of  the  closed-loop  system  is  carried 
out  using  the  reguia tor  control  law  designed  and  the  following  initial 
conditions,  (0(o)  =  0.1  and  <$>(o)  =  0.2]; 


x( 0  )T  --  [0  0  0  .1  0  0  0  .2] 


(62) 


The  responses  of  the  aircraft  to  the  initial  condition  of  Equation  (62) 
after  the  failure  of  the,  right  aileron,  for  the  original  system  are 
shown  in  Figure  12.  In  general,  the  aircraft  responses  to  the 


multivariable  control  law  are  stcble,  settling  time  is  faster,  and 
oscillations  of  the  short  period,  dutch  roll,  and  the  spiral  divergence 
modes  were- removed  from  the  respective  outputs. 

The  eigenvalue  spectrum  of  this  regulator  design  spans  the  Z-plane 
from  0.001  to  0.99.  Even  though  the  fast  eigenvalues  responses  have 
settled,  the  slower  eigenvalues  responses  still  drive  the  system  with 
small  piecewise-constant  control  input  changes  between  the  sampling 
period.  The  magnitude  of  the  control  inputs  did  not  overdrive  the 
system  after  the  initial  sampling  period,  xnus,  a  relatively  smooth 
transition  between  the  sampling  periods  from  the  initial  conditions  to 
the  settling  time  of  shout  4  sec  for  the  pitch  angle  Figure  12a  and  the 
pitch  rate  Figure  12b.  The  fast  settling  time  of  the  roll  angle  in 
Figure  12c  ind  cates  that  the  initial  control  inputs  for  some  of  the 
control  surfaces  were  large.  Consequently,  4>  was  forced  to  approxi¬ 
mately  zero  in  less  than  half  a  sampling  period.  The  left  aileron  and 
spoilers  did  exceed  their  physical  limits.  Also,  the  spoilers  had  a 
negative  and  positive  deflection;  this  is  impossible  since  the  spoilers 
only  have  a  negative  deflection  (CCW)  on  the  A-7D.  The  simulation 
routine  in  CESA  is  for  a  linear  system  and  control  surface  limitations 
are  not  recognized.  The  roll  rate  response  of  Figure  12d  is  well 
within  the  maximum  roll  rate  of  200  deg/sec  (Ref  l).  A  small  initial 
sideslip  angle  (0.0017  rad),  Figure  12e,  and  yaw  rate  (-.012  rad), 
Figure  12f,  does  exist,  but  is  negligible. 

The  block  diagram  for  the  internally  balanced  state-variable 
feedback  controller  is  shown  in  Figure  11.  The  transformed  state 
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Figure  12c.  Roll  Angle,  Initial  Conditions  .2  Rad  vs. 


Figure  12d.  Roll  Rate  Angle  vs. 
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equation  for  the  balanced  system  is,  as  follows: 


x'(t)  =(T":iAT)x'(t)  +  (T-1B)u(t)  (63) 

y(t)  =  CT  X'(t)  (64) 

and  the  control  law  in  the  discrete- time  domain  is 

u(kt)  =  K'x'(kt)  (65) 

The  actual  magnitude  of  the  control  input  is  required  in  order  to 
determine  if  the  designed  control  law  exceeds  the  physical  limits  of 
the  actuators.  Therefore,  since  x'(kt)  =  T  1  x(kt), 

u(kt)  «  K'T-1  x(kt)  =  Kx(kt)  (66) 

Also,  the  initial  condition  inputs  for  the  balanced  system  must  be 
transformed  into  the  new  coordinate  system  by 

x'(o)  =  T”1x(  o )  (67) 

Proceeding  with  the  regulator  design  for  the  internally  balanced 
system  using  Equation  (A-25  through  A-27),  the  closed-loop  eigenvalue 
spectrum  is 

o(F  +GqK)  *  {0.005,  0.007,  0.009,  0.01,  0.03,  0.05,  0.07,  0.09}  (68) 

Assigning  the  eigenvector  for  each  of  Equation  (59),  the  matrix  X' 
of  eigenvectors  is 
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Q.2264E-01 

-0.7282E-01 

-0.3186E-02 

0.3118 

-0.2095E-02 

-0.1041 

-0.9799E-02 

0.1621 

-0.1516 

-0.1033E-01 

-0.3H2E-02 

-0.1016E-01 

0.2754 

0.1845E-U1 

0.1480E-01 

0.2166E-01 

0. 

1.000 

0. 

0. 

0.7869 

0.9943E-02 

-0.3865E-01 

0.1089 

1.000 

0. 

0. 

0. 

0. 

0. 

0. 

1.000 

-0.1738E-03 

0. 3050 

0.3684E-03 

-0.2891E-02* 

-0.9793E-03 

0.1577 

0.9896E-03 

-0.1042E-01 

-0.2651E-0 

-0.1012E-01 

0.9387E-03 

-0.4218E-02 

-0. 5001E-03 

0.2096E-01 

-0.3472E-02 

0.1800E-01 

0. 

0. 

0. 

0. 

0.1549E-02 

0.1077 

0.7428E-02 

-0.3664E-01 

0. 

0. 

0. 

0. 

0. 

1.000 

0. 

0. 

w* 

the  associated 

Q  matrix  is 

"o. 8855 

0. 3501E-01 

-1.241 

-0.2361E-01 

0. 

0. 

1.000 

0. 

0.1126 

0.4050 

-0.2651 

0.9257E-01 

0. 3665 

0.1383 

-0.9529 

-0.2482 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

-0.4475E-01 

-0.U68E-01 

-0.7149E-02 

-1.238  “ 

0. 

0. 

0. 

1.000 

0.4621E-02 

0.9389E-01 

-0.1168E-01 

-0.2644 

-0, 8522E-01 

-0.1959 

0.6394E-01 

-0.9519 

0. 

0. 

1.00 

0. 

1.000 

0. 

0. 

0. 

76 


The  feedback  matrix  of  Equation  (65) 


-131.8 

213.0 

-403.8 

-79.73 

93.20 

-157.2 

376.8 

85.95 

-28.07 

48.09 

-131.6 

-30.82 

-93.77 

129.5 

-64.27 

7.069 

402.0 

-725.2 

2233. 

495.0 

_249.8 

-281.3 

-1757. 

-650.0 

7.978 

-9.034 

-27.82 

5.146* 

-5.473 

8.870 

24.02 

-2.573 

1.960 

-4.226 

-7.283 

.8373 

5.601 

7.432 

-14.77 

6.376 

-22.72 

160.5 

65.12 

-13.27 

-25.46 

-56.75 

-48.79 

-29.86_ 

The  average  gain  for  Equation  (71)  is  199.  Using  the  inverse  of  the 
transformation  matrix  T,  the  matrix  K  given  in  Equation  (66)  is 
obtained  wher«  T"1  is 
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3.1310E-5 

-1.8741E-2 

3.1528E-3 

2.2038E-2 

-1.5149E-5 

1.6794E-2 

-1.9646E-3 

-1 . 8162E-2 

-2.5606E-3 

4.3495E-1 

-4.6C60E-3 

-4.5933E-1 

5.6238E-4 

-7.7507E-1 

1 . 1906E-1 

9.1296E-1 

-3.4657E-5 

1.8281E-2 

-2.9889E-3 

-2.0732E-2 

-1 . 6022E-4 

-5.1974E-4 

2.9203E-1 

-7.4909E-2 

2.8776E-4 

-1.0776 

9.3580E-2 

1.8823E-1 

J-1.4248E-5 

-4.0350E-2 

-3.6914E-2 

-1.4945E-2 

-7.8563E-2 

6.2237E-2 

4.6185E-1 

1.0295 

-2.4268E 

3.4092E-2 

1.S894E-1 

-1.8377 

-7.9595E-3 

4 . 1505E-4 

-3.6800E-3 

1.7059E-2 

-3.8466E-2 

-9.9916E-5 

7.4374E-3 

-8.2169E-2 

6.4536E-1 

-3.4859E-2 

6.8651E-1 

9.1750E-2 

1.2719E-1 

i . 4930E-2 

2.9657E-2 

3.3226E-3 

-8.3743E-2 

-9.2893E-3 

-1.6389E-2 

-S.4078E-3 

1.3260 

1.8622E-1 

3.6654E-1 

-2.2451E-2 

TlUs  results  in  the  actual  feedback  matrix  K  as  illustrated  in  Equation 
(66).  Then  the  actual  controls  can  be  analyzed  to  ensure  that  they  do 
not  exceed  the  physical  limits  of  the  actuators. 


0.9734 

-77.87 

-0.9011 

67.00 

0.3159 

-24.20 

0.1591 

-13.74 

-5.396 

497.7 

4.133 

«*■» 

-216.5 

-487.1 

-0.2977 

360.0 

0.2104 

-110.0 

0.5744E-01 

-292.4 

-0.2085 

1674. 

1.287 

642.8 

0.2244 

2.817 

101.1 

-1.570 

-85.69 

0.2630 

29.75 

3.788 

28.01 

12.30 

-550.8 

-15.08 

220.3 

-42.18 

-520.4 

7.765 

377.8 

-1.053 

-.112.7 

-11.91 

-334.7 

27.34 

1708. 

34.55 

823.1 
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Using  the  K'  feedback  matrix,  the  closed-loop  plant  matrix  has  the 
specified  eigenvalues  of  Equation  (68); 


-2.118 

0.4608 

0.2956E-01 

0.1507 

-1.094 

0.2927 

-0.9719E-01 

-0.1276 

0.8791E-01 

-0.192 5 E-01 

0.2463E-01 

-0.8952E-01 

[F'+G'K']  = 

-0.2155 

0.3584E-01 

0.6555E-01 

0.4149 

-0.5578E-12 

0.4317E-12 

-0.1007E-11 

0.1473E-11 

-0.6021 

0.1319 

-0.6959E-01 

-0.7137 

0.9511E-12 

-0.1619E-11 

0.3949E-11 

0.70172-13 

-6.994 

1.482 

0.1815 

0.7074 

-0.1102 

0.1138 

-0.7759E-01 

0.5732 

-0.4885E-01 

-0.2275E-01 

0.6369E-01 

0.2996 

0.6697E-02 

-0.3791E-01 

0.5543E-01 

-0.1807E-01 

-0.2005E-01 

0.1588 

-0.2230 

0.3598E-01 

0.7000E-02 

0.3382E-12 

-0.9993E-12 

0.1066E-13 

-0.1553E-01 

-0.2425 

0.3947 

0.2086 

-0.7200E-13 

-0.2984E-12 

0. 5100E-02 

-0.2665E-14 

-0.3699 

0.4572 

-0.3656 

1.887 

Before  the  continuous-time  simulation  is  carried  out,  the  initial 


conditions  of  Equation  (62)  are  transformed  by  Equation  (67)  to 


x'(o) 


0.20810 

-0.36936 

-0.04218 

0.07486 

0.01628 

-6.826E-3 

0.817141 

-5.9348E-3 


(75) 


The  continuous  time  ^lots  of  the  balanced  system  are  shown  in 
Figure  13.  The  eigenvalue  spectrum,  Equation  (68),  for  this  regulator 
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design  was  located  close  to  the  origin  in  the  Z-plane.  These  eigen¬ 
values  have  a  faster  transient  time  than  the  previous  design.  Thus, 
the  control  law  realised  drives  the  response  of  the  system  harder 
between  the  sampling  periods  to  a  settling  time  of  about  4  sec  causing 
what  appears  to  be  the  short  period  or  dutch  roll  oscillation.  However, 
that  is  not  the  case.  The  increased  magnitude  of  the  piecewise-constant 
control  inputs  are  overdriving  the  system,  producing  the  oscillations 
and  not  the  aircraft  modes.  An  important  result  in  all  the  balanced 
responses  is  that  the  initial  conditions  are  approximately  the  same  as 
for  the  original  system.  Thus,  the  transfoi'mation  of  Equation  (67)  on 
the  initial  conditions  and  the  internally  balanced  matrices  do  represent 
the  original  system.  The  responses  differ  from  the  original  system 
because  of  a  different  eigenvalue  spectrum  and  eigenvectors.  In 
general,  disregarding  the  overshoots  produced  by  the  increased  constant 
piecewise-constant  inputs,  the  responses,  Figure  13,  do  show  an  improve¬ 
ment  in  the  settling  time  of  the  system. 

The  internally  balanced  matrices  of  Equation  (A-25  through  27)  were 
reduced  to  seven  states  as  illustrated  in  Equation  (43).  Using  the 
following  closed-loop  eigenvalue  spectrum: 

°(F,R+0RKr)  =  {0.005,  0.007,  0.009,  0.01,  0.03,  0.05,  0.07}  (76) 

The  eigenvectors  of  Equation  (57)  were  assigned  for  each  X^  of  Equation 
(76)  and  resulted  in  the  following  eigenvector  matrix  X' 
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0.2293 

0.5702E-02 

-0.5106 

-0.5323 

0.1166 

-0.5313E-01 

-0.3096 

-0.3148 

il 

\  cs 

X 

-0.1502 

-0.1484E-01 

-0.7933E-02 

-0.5111E-02 

0.2752 

0.1838E-01 

0.1839E-01 

0.3365E-02 

0. 

1.000 

0. 

0. 

0.8248 

0.2442E-01 

-0.1227 

-0.8793E-01 

1.000 

0. 

0. 

0. 

-0.4551E-01 

-0.5428 

0.3549E-01 

-0.2850E-01 

-0.3238 

0.22C1E-01 

-0.7559E-03 

-0.6274E-02 

0.1376E-02 

-0.1669E-03 

0.6174E-02 

-0.3136E-02 

0. 

0. 

0. 

-0. 5911K-02 

-0.8523E-01 

0.1282E-01 

0. 

0. 

0. 

and  the  associated  ma 

trix  u  of  the  Equation  (59)  is 

— 0.778J 

-0.5428E-02 

-0.9561 

0.2996 

0. 

0. 

1.000 

0. 

V 

0.1013 

0.4008 

-0.2206 

0 . 4709E-01 

0. 

0. 

0. 

1.000 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

-0.1872^-01 

0.3042 

-0.2672E-01 

0. 

0. 

0. 

0.9154E-02 

9.5332E-01 

-0. 14753-01 

0. 

1.000 

0. 

0. 

0. 

1.000 

1.000 

0. 

0. 
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The  matrix  of  Equation  (65)  is 


-106.5 

176.1 

-294.3 

23.13 

78.12 

-131.7 

299.4 

12.99 

-23.28 

39.44 

-105.1 

-5.832 

-68.78 

105.8 

-0.4973 

66.33 

329.0 

-588.0 

1S10. 

97.25 

144.7 

-181.6 

-2025. 

-899.4 

5.335 

29.11 

-69.93 

-3.452 

-IS. 26 

53.91 

1.250 

5.071 

-17.52 

4.506 

29.13 

-33.92 

-11.28 

12.40 

221.1 

-21.29 

-147.6 

53.05 

The  olosed-l cop  discrete  plant  matrix  of  the  reduced  balanced 
system  for  the  ass  lined  eigenvalues  is 


-0.7759 
-0.4539 
0.9493E-02 
-0.6  JV3K-01 
|  0.58622-12 
|  0.2704E-01 


1.965 

1.054 

-0.103;E-0l 

0.1065 

-0.9134E-12 

0./:'03E-01 


-6.184 

-3.565 

0.5284E-01 

-0.2890 

0.4164E-11 

-0.2270 


m0. 95322-12 

-0.1353E-11 

0.3867E-1] 

0. LS81 

-2.174 

2.433 

0.1121 

-] .294 

1.462 

0. 1652E--32 

-0.2;80E~01 

0.419:2-01 

0.5358E-03 

0.24 381-01 

-0.73  '42-01 

o.voooE-oa 

0.3864E-13 

-Q.7jf‘-E-12 

0.2 382 E-01 

-0.3303 

0.4811 

-0.5L18E-13 

-0.3197E-12 

0 . 5000E-02 

-5.084 

-3.491 

-0.5843E-01 

0.68670-01 

0.1756E-11 

-0.9062 

0.5089E-12 
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The  initial  condition  vector  x'(o)  of  Equation  (75)  was  reduced  by  re¬ 
moval  of  the  last  element  for  the  continuous -time  simulation.  Figure 
14  is  the  time  response  plot  for  the  reduced  balanced  system  and  is 
compared  to  the  previous  cases  for  improved  or  deteriorated  responses. 

The  regulator  design  for  the  reduced  balanced  system  exhibits 
approximately  the  same  initial  conditions  as  the  two  previous  designs. 
Due  to  the  fast  transient  responses  of  the  eigenvalue  spectrum,  the 
feedback  matrix  produces  a  control  law  that  overdrives  the  system 
during  the  first  two  sampling  periods  only.  The  pitch  angle,  Figure 
14a,  differs  from  Figure  13b  only  little  in  the  initial  response  and 
the  settling  time  is  decreased  from  4.0  sec  to  1.5  sec  since  the  control 
input  induced  oscillation  is  not  present.  For  the  remainder  of  the 
outputs,  the  settling  time  decreased  from  4.0  sec  to  2.5  sec.  By  re¬ 
moving  the  roll  subsidence  state,  which  had  a  low  degree  of  controlla¬ 
bility,  the  aircraft  model  was  not  degraded.  Thus,  the  magnitude  of 
the  control  inputs  was  decreased  sufficiently  to  remove  the  oscillation 
caused  by  the  control  inputs  of  the  balanced  system. 

To  transform  the  reduced  balanced  feedback  matrix  (K£)  in  the  new 
coordinate  system  back  to  the  original  coordinate  system  the  following 
transformation  is  applied.  Firs+  the  reduced  balanced  state  vector  is 
obtained  as  follows, 


x'  (t)  «  fl  .  ,  0  .  .1  x'(t) 

R  l  n-j,  n-j ,  n-j,  jJ 


where  j  is  the  number  of  balanced  states  to  be  removed,  n  is  the 
dimension  of  the  original  plant  matrix,  I  is  an  identity  matrix,  and 
j  is  the  number  of  zero  (0)  vectors.  The  state  equations  for  the  reduced 


balanced  system  are, 


i£(t)  «  A£  x£(t)  +  B'  u(t) 

R 


where  the  matrices  of  Equation  (82)  are  obtained  as  shown  in  Equation 
(43)  and  x'  eR7.  The  eigenvalues  of  the  reduced  model  are  shifted, 
hopefully  only  by  a  small  amount  and  still  retaining  the  characteris¬ 
tics  of  the  model.  The  control  law  for  the  reduced  system  is 

u(t)  =  K£(t)  x'(t)  (83) 

where  is  of  dimension  wx(n-j),  and  w  is  the  column  dimension  of  the 
control  input  matrix.  The  transformation  of  the  control  law  to  the  bal¬ 
anced  system,  using  Equation  (81)  is 

K-J,  n-J,  Vj,  j] 

To  obtain  the  actual  K  matrix  in  the  original  state  coordinate  system, 
x'  =  T-1x  is  substituted  into  Equation  (84)  to  form 

u( t )  s  (yJ*  f I  ,  .  0  .  .1  T'1  |  x( t )  (85) 

L  n-j,  n-j,  n-j,  jj  nxnl  ' 

Therefore, 


u(t)  =  K  x(t) 
wxn 


once  the  transformation  of  the  reduced  balanced  system  control  law  to 
the  original  system  is  completed.  The  time  responses  produced  by  the 
transformed  control  law  can  then  be  compared  to  the  original  time 
responses.  This  allows  one  to  investigate  the  effects  of  using  a 
reduced  feedback  matrix,  transformed  back  into  the  original  state 


i 


coordinate  system  in  the  control  laws. 


Summary 

By  increasing  the  control  input  capabilities  of  the  CESA  arrays 
and  determining  the  degree  of  controllability  of  the  aircraft  states, 
a  multivariable  control  law  was  designed.  The  time  response  plot-  of 
the  original,  and  internally  balanced  systems  proved  that  the  state 
coordinate  system  can  be  changed  without  modifying  ti  e  sj/"tem 
characteristics.  The  control  laws  designed  for  a  sampled  data 
regulator  and  tracker  illustrate  the  capabilities  of  the  entire  eigen- 
structure  assignment  design  method.  The  assignment  of  the  eigenvalue 
spectrum  and  associated  eigenvectors  points  out  their  effect  on  the 
closed-loop  system  as  shown  in  the  time  response  plots  for  the 
different  regulator  cases.  Due  to  the  large  difference  in  the.  eigen¬ 
value  spectra  used  in  the  regulator  designs,  the  system  responses 
exhibit  an  acceptable  settling  time  and  a  decrease  in  the  overshoots. 
The  desirable  feedback  matrix  should  have  an  average  gain  about  1.0, 
due  to  the  eigenvectors  selected,  and  would  constitute  an  acceptable 
control  law. 


00a 


REDUCED  BALANCED  SYSTEM  RIGHT  AILERON  FAILURE 

CONTINUOUS  SIMULATION  SYSTEM  OUTPUT  »  1 


o 


(SECONDS) 


REDUCED  BALANCED  SYSTEM  RIGHT  AILERON  FAILURE 


SECONOS) 


% 


Figure  14f.  Yaw  Rate 


V.  Conclusions  and  Recommendations 

Conclusions 

Reconfiguring  the  control  surfaces  response  to  compensate  for  a 
primary  surface  failure  required  a  comprehensive  aircraft  model. 

Using  the  left  and  right  control  surfaces  as  independently  maneuverable 
inputs  in  the  six  degree-of-freedom  aircraft  equations  of  motion#  an 
accurate  aircraft  model  wa3  developeu  in  Chapter  II.  The  importance 
of  coupling  all  the  control  surface  inputs  to  both  the  longitudinal 
and  lateral  axes  is  illustrated  by  the  time  responses  for  the  uncon¬ 
trolled  original  system  in  Chapter  III.  Therefore,  the  equations  of 
motion  cannot  be  decoupled  when  a  sampled-data  controller  is  designed. 
Inclusion  of  this  coupling  between  the  axes  required  the  derivation  of 
several  new  non-dimensional  control  derivatives.  The  new  dimensional 
control  derivatives  were  then  calculated  for  use  in  the  aircraft 
model.  Using  the  geometrical  properties  of  the  aircraft  and  the 
Digital  Datcora  computer  program,  the  needed  control  derivatives  were 
derived  as  shown  in  Appendix  B. 

A  comprehensive  review  of  entire  eigenstructure  assignment  is 
presented  in  Chapter  III.  The  use  of  this  method  for  designing 
multivariable  control  laws  is  feasible  since  there  is  complete  freedom 
in  assigning  the  eigenvalues  and  the  eigenvectors  must  be  within 
specified  sub  spaces.  With  a  comprehensive  model  now  available,  the 
entire  eigenstructure  assignment  method  was  used  to  assign  both  the 
eigenvalues  and  eigenvectors  to  the  closed-loop  plant  matrix  by 
synthesizing  a  state  variable  feedback  control  la w  for  a  regulator 
and  tracker.  Using  the  interactive  computer  program  CESA,  a  iirect 
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digital  design  of  both  a  sainpled-data  regulator  and  tracker  were 
attempted  in  Chapter  IV.  Considering  the  number  of  control  inputs, 
the  degree  of  controllability  of  the  states  was  established  by  the 
singular  value  decomposition  of  the  H^p  matrix  by  means  of  Moore’s 
algorithm  for  internally  balanced  matrices.  Model  reduction  is  then 
possible  by  removing  the  state(s)  of  low  controllability  as  indicated 
by  the  magnitude  of  the  singular  values. 

With  the  degree  of  controllability  of  the  system  determined,  the 
regulator  control  law  design  should  be  carried  out  first  to  uncover 
any  mathematical  pi'oblems  with  a  new  system.  The  regulator  designs 
for  reconfiguring  the  control  surfaces  of  this  aircraft  model  results 
in  a  aero  steady  state  error.  However,  the  "proper"  selection  of  the 
eigenvectors  would  improve  the  magnitude  of  the  initial  control  inputs 
so  that  the  physical  limits  of  the  actuator  are  not  exceeded.  Were  it 
not  for  the  time  limit  at  this  point  in  the  thesis,  it  would  be  possible 
to  reimplement  the  tracker  design  for  a  reeonfigurable  control  law. 

Plots  of  the  regulator  designs  are  included  in  this  thesis. 

Recommendations 

The  aircraft  model  developed  in  this  thesis  and  the  attempt  to 
design  a  sampled-data  controller  by  the  entire  eigenstructure  assign¬ 
ment  method  paves  the  way  for  further  studies  in  the  area  of  recon- 
figurable  control  laws .  Areas  that  need  further  investigation  and 
development  are: 

1.  During  the  design  of  an  aircraft,  more  extensive  wind  tunnel 
testing  is  necessary  to  provide  the  data  for  the  additional  control 


93 


derivatives  in  order  to  obtain  a  more  accurate  aircraft  model  driven  by 
individual  control  surfaces. 

2.  The  development  of  the  criteria  for  the  selection  of  the 
eigenvalues  and  eigenvectors  assigned  to  the  closed-loop  plant  matrix 
to  insure  the  design  of  a  realizable  multivariable  control  law.  This 
includes  the  further  study  of  the  entire  eigenstructure  assignment 
method  to  modify  the  control  surface  input  to  compensate  for  a  parti¬ 
cular  control  surface  failure. 

3.  The  development  of  a  fault  and  isolation  detection  system  that 
monitors  the  control  surface  responses.  This  would  determine  which 
control  surface  is  not  following  the  command  inputs.  This  may  be 
accomplished  by  means  of  a  Kalman  filter. 

4.  The  investigation  of  a  digital  design  method  that  uses  the 
plant  end  the  measurement  (sensor)  outputs  as  feedback  which  are 
compared  with  the  input  commands  for  the  design  of  an  error-actuated 
controller  as  it  uould  apply  to  reconfigurable  tracker  control  laws. 

5.  The  modification  of  the  interactive  program  CESA  to  facilitate 
the  transformation  of  the  reduced  feedback  balanced  matrix  to  the 
original  state  coordinate  system.  Then  apply  the  transformed  state 
feedback  to  the  original  system  and  compare  its  response  with  that  of 
the  reduced  system. 
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APPENDIX  A 


Introduction 

This  appendix  presents  a  model  of  the  A-7D  using  the  linearized 
aircraft  equations  of  motion.  Coupling  between  the  longitudinal  and 
the  lateral-directional  axes  is  through  the  control  surfaces  described 
earlier.  The  linearized  differential  equations  of  motions  are  then 
used  to  derive  the  continuous  state  space  model. 

Equations  of  Motion 

The  lateral-directional  and  longitudinal  equations  are  developed 
with  coupling  between  the  axes  through  the  non-traditional  control 
inputs  (Ref  4).  These  equations  assume: 

1.  X,  Y,  and  Z  axes  are  in  the  plane  of  symmetry  and  the  origin 
of  the  axes  is  at  the  center  of  gravity  of  the  aircraft. 

2.  The  mass  of  the  aircraft  is  constant . 

3.  The  aircraft  is  a  rigid  body. 

4.  The  earth  is  an  inertial  reference. 

5.  The  perturbations  from  equilibrium  are  small. 

6.  The  flow  is  quasisteady  (Ref  2).  The  linearized  longitudinal 
equations  of  motion  are  as  follows: 

u  «  -gecose,  +  Xu  +  X  a  +  X»a  +  X.  6,,  ♦  X.  «  +  X,  6„ 

*  u  a  a  6^  it  r  a 

w  =  Uiq-gesinGi  ♦  Z  u  +  Z  a  +  Z*a  +  Z  a  +  Z.  6..  +  Zx  6 


tho  dimensional  stability  derivatives  for  the  aircraft  model  are 
given  in  Table  I.  Excluded  are  the  dimensional  control  derivatives 
which  are  derived  in  Appendix  B.  All  of  the  derivatives  relating  to 
angles  are  per  radian  measure  (Eef  21).  The  continuous  state  space  of 

the  model  is  now  formed. 


Continuous  State  Space 

The  aircraft  model  has  eight  state  variables  and  the  control 
inputs  are  defined  as  pairs  as  shown  below.  When  used  as  independent 
control  surfaces,  «u  becomes  the  5^  (left)  and  (right)  horizontal 

stabilizer.  The  remaining  sets  of  control  surfaces,  except  the  rudder, 
are  split  into  the  left  and  right  control  surface  when  the  multi- 
variable  control  laws  are  designed.  The  states  are: 

1.  u  Perturbation  Forward  Velocity 

2.  a  Perturbation  Angle  of  Attack 


3,  q  Perturbation  Pitch  Pate 

4,  8  Perturbation  Pitch  Angle 

5,  8  Perturbation  Sideslip 

6,  p  Perturbation  Roll  Rate 

7,  r  perturbation  Yaw  Rate 
S,  $  Perturbation  Roll  Angle 

The  control  inputs  are: 

1.  <5.  Horizontal  Stabilizer  Deflection 
it 

o  “  Rudder  Deflection 

r 

3.  6  Aileron  Deflection 

a 

4.  5  Spoiler  Deflection 

s 

5.  6  Flap  Deflection 
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The  dimensional  stability  derivatives  for  the  aircraft  model  are 
given  in  Table  I.  Excluded  are  the  dimensional  control  derivatives 
which  are  derived  in  Appendix  D.  All  of  the  derivatives  relating  to 
angles  are  per  radian  measure  (Ref  21),  The  continuous  state  space  of 
the  model  is  now  formed. 

Continuous  State  Space 

The  aircraft  model  has  eight  state  variables  and  the  control 

inputs  are  defined  as  pairs  as  shown  below.  When  used  as  independent 

control  surfaces,  6^  becomes  the  5h  (left)  and  5^  (right)  horizontal 

1  r 

stabilizer.  The  remaining  sets  of  control  surfaces,  except  the  rudder, 
are  split  into  the  left  and  right  control  surface  when  the  multi- 
variable  control  laws  are  designed.  The  states  are: 

1.  u  Pertui'bation  Forward  Velocity 

2.  a  Perturbation  Angle  of  Attack 

3.  q  Pertui'bation  Pitch  Rate 

4.  8  Perturbation  Pitch  Angle 

5.  B  Perturbation  Sideslip 

6.  p  Pertui'bation  Roll  Rate 

7.  r  Perturbation  Yaw  Rate 

8.  <J>  Perturbation  Roll  Angle 
The  control  inputs  are: 

1.  6..  Horizontal  Stabilizer  Deflection 

it 

2.  C  Rudder  Deflection 

r 

3.  <5  Aileron  Deflection 

a 

4.  5  Suoiler  Deflection 

s 

5.  Flap  Deflection 
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TABLE  I 


Cruise  Configuration 
Dimensional  Stability  Derivatives 


X 

-0,00829  * 

u 

X 

5.47751 

0 

g 

32.2  ft/sec2 

z 

-0.11324 

u 

z 

-632.63786  I 

0 

z 

0.0 

* 

M 

0.00036 

u  1 

M 

-8.15547 

a 

M* 

-0.11565 

M 

q 

-0.59329 

*8 

-102.69611 

Y 

0.56427 

P 

Y 

1.35163  i 

r 

L8 

-25.73501 

L 

-3.00146 

P 

L 

0.90217  i 

r 

3.99375 

N 

-0.0067 

P 

N 

-0.5096 

r 

xunits  rad 


Flight  conditions  for  Table  I  are  given  in  Appendix  B,  Table  I. 


The  continuous  state  space  equations  of  motion  are  derived  from 
Equations  (A-l  -  A-6).  The  state  equations  are  then  developed  for  use 
with  the  CESA  interactive  computer  program  CESA,  so  that  a  discrete¬ 
time  state  variable  feedback  control  law  can  be  designed. 

To  derive  the  state  space  equations,  one  must  solve  for  the  time 
derivatives  of  u,  a,  q,  0,  0,  r,  p,  and  <J>.  This  derivation  assumes  Xct, 
Za,  and  Sq  are  zero  (Ref  21).  Also,  for  a  steady  level  flight  condi¬ 
tion,  Uj  =  constant,  Vj  =  0,  Wj  =  0,  using  the  stability  axis,  0*  =  0, 

=  0,  and  Pj  -  Qi  =  Ri  -0. 

Equation  (A-l)  solved  for  u: 


u  =  -g9  +  Xu  +  X  o  +  X.  <S. . 

u  °  sn  u 


X,  6  +  x.  6  +  X.  6  +  X.  6« 
°r  r  6  a  6  s  6.  f 
a  s  f 


(A-9) 


Equation  (A-2)  is  solved  for  a: 

Letting  w  =  ua  implies  w  =  U^a  +  u  04  where  uaj  -  0.  Since  aj at 
equilibrium  is  aero,  therefore,  w  -  Uj  a. .Substituting  into  Equation 
(A-2)  yields 


Uict  =  Uiq  +  Z  u  +  Z  a  +  Z  6,  +  Z.  6 

*  u  u  5.,  it  6  r 

it  r 


+  Z.  6n  +  Z.  6o  *  Z.  6- 
oa  6  s  6„  f 

U  S  1 


( A-10) 
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Equation  (A-3)  is  solved  for  q. 


4  =  V  +  V  +  +  V  +  +  \5r 


+  Va  +  \*b  +  M6  *f 


(A-12) 


To  express  q  in  terms  of  Ujci,  q  Equation  (A-ll)  is  substituted 
into  Equation  (A-12)  for  u,  then  regrouping  the  terms} 

q  ■  M^u  +  Mqo  +  Ma(q  +  Zu  U  +  Za<x)  ♦  Mqq 

Uj  Uj 


♦A1 Vs  T +  \Vf 


Y^fiit *  m64N  fiit  y  4  \yr  +f  \ 4  $a 


(A-13) 


The  final  equation  needed  in  the  longitudinal  state  space  representation 


0  *  q 


(A-14) 


Equation  (A- 5)  is  solved  for  B: 

Let  V  s  Uj  I,  which  is  substituted  into  Equation  (A-4) 


8  *  i  ^ee  *  et  *  ♦/  Vi 


*lE 

u,l 


v*  ‘ v*  *  v* 


•  V  •  v>] 


(A-15) 
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Equation  (A-6)  is  solved  for  r: 

The  EL  =  p  -I  r  and  to  simplify  use,  a  new  set  of  primed  stability 
I 

xx 

derivatives  (Ref  6)  are  used: 


I 

xz 


L,  +  I  N,  . 
i  xx  i  and 


1  -  V 


WV 


H.'  *  H  *  —  L. 
i  i  I  i 
zz 


1  - 1. 


XZ 


I  (I  ) 

XX  zz 


(A-16) 


where  i  represents  8,  p,  r,  6it'  5a'  ^r'  5s'  and  5f’ 

#  ■  v  *  *  v  *  *  v* 

*  L5  Sa  *  lS  Sb  *  LLSf  (A-17) 

a  s  1 

Equation  (A-7)  is  solved  for  p: 

The  EN  =  r  -I  p  and  simplify  Equation  (A-6)  as  above. 

XZ 

I 

zz 


p  =  N'8  ♦  M'p  +  M'r  +  H'  6.,  +  N'  « 
y  3  P'r  6it  it  5p  r 


+  m;  6  +  n:  «  *  h;  s, 

6a  a  5S  3  6f  «• 


(A-18) 


The  last  equation  for  the  lateral  axis  equation  is: 
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With  a  5  failure  the  B matrix  is: 
ar 


-16.43 


!  -7.955 


-16.43 


-0.6730E-1  -0.6730E-1 


-7.955 


-7.506 


-0.4788E-1 

1.860 


-8.040 

0.6014 


8.040 

-0.6014 


0.4530E-1  -0.5816E-2 


5.967 

-5.198 


17.19 

0.2985E-1 


-0.2013 


-1.260 


-0.2013 


-9.3423 


-0.9707E-2  -0.9707E-2  -0.9351E-1 


-1.0488 


-9.3423  «hl 
-0.9351E-1  6^ 
-1.0488  5r 


0.1424E-2 


-1.099 


-0. 5335E-1 


-0.1424E-2 

1.099 


0.5335E-1 


-5.3212 


0.1571 


5.3212 

-0.1571 


#See  page  116. 


(A-23) 


where  the  failed  right  aileron  control  input  is  missing.  The  sign  con¬ 
vention  for  the  control  surface  is:  positive  6^,  an<*  5f  *s  trailing 
edge  down;  positive  6r  is  trailing  edge  left;  positive  <5S  is  spoiler  up. 
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For  the  regulator  design  the  left  and  right  flap  inputs  are  dropped  to 
resolve  the  controllability  and  eigenvector  problems  associated  with 
the  tracker  design. 

The  tracker  output  equation  is: 


0  -1  010000 


Cx  = 


0  0001000 
0  0000010 


0  0000001 


q 

9 

6 

p 


r 


( A-24 ) 


and  the  desired  command  input  for  the  tracker  is: 


1 


°  0-a  3  Y 


v(t) 


0  3  8 

0  ~  r 


(A-25) 


where  Y  =  0-a  is  the  flight  path  angle.  The  regulator  output  matrix  C 
is: 


00100000 

•  - 

u 

00010000 

a 

q 

00001000 

0 

00000100 

8 

00000010 

p 

r 

00000001 

A 

(A-26) 
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IWKa*. 


aSWSSggBB 


The  following  equations  are  for  the  internally  balanced  matrices 
of  the  system  generated  by  MBS)  (Ref  13).  The  plant  matrix  is: 


-0.1695 

-0.2289E-2 

-0.1297E-1 

0.1447E-1 

-0.2146 

-0.7590E-4 

-0.6469E-2 

0.1587E-1 

0.2143E-1 

0.9030E-2 

0.5677E-1 

0.1149 

-0.3819E-1 

-0.2011E-1 

-0.7662E-1 

-0.8726E-1 

-1.663 

-1.004 

-0.2189E-1 

0.5202E-1 

-0.6214E-1 

-0.2490E-1 

0.2280 

-0.5543 

-0. 1076E-1 

-0.2002E-2 

0.3949E-1 

-0.2044 

-0.2895 

-0.1751 

0.1538E-1 

-0.3358E-1 

1.718 

-0.5732E-1 

0.4076E-1 

-0.4794 

1.025 

-0.3078E-1 

0.1075E-1 

-0.2376 

0.2811E-1 

-0.3121E-1 

0.7517E-2 

-0.8604E-2 

-0.4492E-1 

0.1865 

0.2124 

0.7503E-2 

-0.7259 

0.1813 

-0.1018 

1.360 

-0.2004E-1 

-1.450 

2.701 

0.4078E-1 

0.3793E-1 

-3.051 

-0.193 

0.3932 

-0.2407 

-0.3918 

-0.1836 

-2.784 

t 
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The  B'  input  matrix  is: 


-0.2470 

0.1983 

-0.1457 

0.1752 

0.3737 

0.3848 

-0.8989 

-0.9095 

0.7163 

-0.6700 

-2.423 

-2.218 

-0.6118 

-0.7415 

-0.9802 

1.573 

0.9265E-1 

1.091 

0.5125E-1 

0.6013 

0.8547E-2 

-0.9944E-1 

-0.1681E-1 

0.2573 

-0.2134E-2 

-0.5889 

-0.4078E-1 

0.3001 

-0.1970E-1 

0.6600E-1 

0.2301 

3.137 

-2.033 

-0.8386 

0.4047E-1 

-0.4522E-1 

-3.747 

-0.5811E-1 

0.2381E-1 

0.7340 


-0.9363E-1 

-0.5074E-1 

0.80Q5E-2 

-0.1749E-1 

0.3069E-2 

-0.7639E-1 


0.2221E-2 


-0.2145 
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The  C'  matrix  is: 
o 


■> 

0.2767E-1 

0.1357E-1 

-0.1832 

0.2226 

-0.4223E-1 

-0.3246E-1 

0.2157 

1.366 

-0.4086 

-0.2212 

-0.7995E-2 

0.1546E-1 

2.183 

0.9912 

0.2776E-3 

0.3343E-1 

0.4128 

0.2848 

-0.7939E-2 

0.1697E-1 

0.6222 

«■» 

-0.2040 

0.1058E-1 

-0.2820E-1 

-0.1859E-1 

3.339 

-0.2253 

-0.2936 

-0.2395E-1 

-0.2415 

-0.8955 

-0.3772E-1 

0.2191 

0.2959E-1 

-0.1504E-1 

0.2150 

-3.747 

0.4795 

-0.2491 

3.707 

1.093 

0.2318E-2 

-0.101E-3 

0J423E-1 

-0.2464 

-0.2907E-1 

0.1439E-1 

-0.2137 

(A-29) 


After  final  review  of  the  B  matrix,  it  was  determined  that  the 
signs  of  the  lateral-directional  coefficients  5hj.»  5f^>  5frJ 

and  the  longitudinal  coefficients  (pitching  moment)  and  6S^, 

<$s ^  (lift)J  required  a  sign  change.  Also  the  6a^  (pitching  moment) 
coefficient  was  corrected.  The  time  response  plots  for  this  thesi.F 
were  obtained  using  the  B  matrix  of  Equation  (A-23).  The  following  B 
matrix  should  be  used  for  future  studies  of  reconfigurable  multivariable 
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control  laws 


Jft 

i 

i> 


I 


* 


C 


-16.43 

-16.43 

-7.506 

0 

-0.6730E-1 

-0.6730E-1 

0 

-0.4788E-1 

-7.955 

-7.955 

0 

-0.3944 

0 

0 

0 

0 

0 

0 

0.4530E-1 

-0.5816E-2 

8.040 

-8.040 

5.967 

17.19 

-0.6014 

0.6014 

-5.193 

0.2985E-1 

0 

0 

0 

0 

-1.260 

-1.260 

-9.3423 

-9.3423 

0.9707E-2 

0.9707E-2 

-0.9351E-1 

-0.9351E-1 

-0.2013 

-0.2013 

-1.0488 

-1 . 0433 

0 

0 

0 

0 

0.1424S-2 

-0.1424E-2 

0 

0 

-1.099 

1.099 

5.3212 

-5.3212 

-0.5335E-1 

0.5335E-1 

-0.1571 

0.1571 

0 

0 

0 

0 

( A- 30 ) 
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APPENDIX  B 


Introduction 

At  this  point,  one  should  remember  that  the  control  surface  inputs 

will  be  separated  into  independent  variable  left  and  right  control 

surfaces.  When  this  is  done,  caution  must  be  exercised  in  assigning 

the  correct  sign  to  the  left  and  right  oontro*  derivative.  For 

example,  the  derivative  L'  ,  the  change  of  rolling  moment  due  to 

°a 

aileron  deflection,  is  -34.4  for  both  surfaces.  This  indicates  that 
the  right  aileron  is  down  and  the  left  aileron  is  up  causing  a  rolling 
moment  to  the  left.  When  the  aileron  surfaces  are  treated  as  inde¬ 
pendent  controls,  the  rolling  moment  coefficient  due  to  the  right 

aileron  (L'  )  has  a  value  of  -17.2  and  the  left  aileron  17.2.  This 

°ar 

a 8 surges  the  sign  convention  is  trailing  edge  down  and  represents  a 
poeitive  deflection. 

This  appendix  presents  the  equations  used  to  derive  the  dimensional 
control  derivatives  and  the  derive tion  of  the  non-dimensional  control 
derivatives  for  non-traditional  control  surfaces  used  as  input  in  the 
aircraft  equations  of  motion. 

Dimensional  Control  Derivatives  Equations 

The  equations  for  the  longitudinal  dimensional  control  derivatives 
are,  as  follows: 

Let  X  =  qS/m  then; 


(B-l) 


Where  5  =  <5..,  5  ,  6  ,  6  ,  <5.,  which  remains  the  same  tor  the  remainder 
1 1  r  a  s  x 
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of  the  equations  illustrated. 

For  the  Z  force  inputs,  let  Z  =  qS/m  then, 


Z  „  * 


ZCL. 


and  for  the  M  torque  inputs,  let  M  3  qSF/I , 


L  =  MC  . 
®  ra$ 


(B-2) 


(B-3) 


Continuing  with  the  lateral-directional  control  derivatives,  let 
Y  3  qS/m 


(B-4) 


and  for  the  L  and  N  torque  inputs,  let  L  3 


qSb/I  „  then, 
zz 


L 


5 


(B-5) 


and 


(B-6) 


The  data  for  use  with  the  above  equations  is  given  in  Table  I,  (Ref  21). 
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TABLE  I 


Aircraft  Data  Cruise  Configuration 


ALTITUDE 

15,000  ft. 

MACH 

.6 

WEIGHT 

25,338  lbs. 

c.g. 

28.71JS  of  mgc 

q 

300.88  lbs/ft? 

8 

375  ft? 

b 

38.73  ft. 

c 

10.84  ft. 

15,365  Slug  ft? 

*zz 

79,005  Slug  ft? 

V 

69,528  Slug  ft? 

rxz 

-1,664  Slug  ft? 

Now  the  method  mentioned  in  Chapter  II  is  used  to  obtain  the  new  non- 
dimensional  control  derivatives. 

Derivation  of  Non-Dimensional  Control  Derivatives 

The  control  derivatives  of  Table  II  are  for  the  aircraft  at  a 
cruise  configuration  of  0.6  Mach,  angle  of  attack  at  4.2°  and  flying  at 
15K  feet  from  the  following  three  references:  digital  DATCGM  (Ref  14), 
A-7  Aerodynamic  Data  (Ref  11),  and  from  (Ref  21).  The  digital  DATCOM 
program  was  used  to  find  those  non-dimensional  control  derivatives  not 
available  from  the  other  reference  sources.  Therefore,  only  those 
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TABLE  II 


Non-Dimensional  Control  Derivatives* 


LONGITUDINAL 

LATERAL-DIRECTIONAL 

cx 

0.1146 

s. 

0.0 

s 

0.29796 

CU 

«r 

-0.02827 

°m5h 

nr 

-0.45276 

nr 

-0.01087 

°*r 

0.05233 

°y« 

ar 

0.02505 

S 

0.0 

X 

-0.06045 

Cmj 

°r 

0.0 

X 

0.0071 

X 

0.0 

cy«r 

0.20055 

X 

0.2120 

V 

0.01902 

X 

-0.0229 

s 

-0.0917 

%s 

sr 

8.7859E-3 

Cy<5s 

sr 

-6.303E-3 

°L«s 

sr 

-0.04297 

Cl«sp 

3.8636E-3 

°sr 

-0.01146 

X 

1.9100E-3 

°Psf 

lr 

0.06515 

X 

0.0 

%f 

lr 

0.41399 

s 

lr 

-0.01865 

lT 

-0.06024 

Sr 

ir 

8.0879E-4 

* units  rad  ,  right  control  surface,  except  rudder 


ri 


1 


f  *4 
r 

'2 


i 

p* 


:v> 


control  derivatives  found  with  digital  DATCOM  will  be  addressed  here. 

The  following  derivatives  were  derived  using  digital  DATCOM:  Cjj  , 

5h 

Cp  »  Cj)  ,  Ct  ,  C-  ,  Ci  ,  C «  and  Ci  since  this  data  is  not 

op  6f*  of*  4  of*  of*  Of 

usually  available  from  wind  tunnel  testing.  To  establir’  confidence 

in  this  method,  a  comparison  was  done  on  between  digital  DATCOM 

6s 

and  the  A-7  aerodynamic  data,  which  is  shown  later  in  this  Appendix. 
For  the  longitudinal  mode,  digital  DATCOM  provided  tha  data  for  the 
following  non-dimensional  control  derivatives:  Cp  ,  Cp  ,  Cp  , 

Cj^  ,  and  Cj^  .  The  drag  for  the  horizontal  stabilizer  is  usually 

considered  zero.  But  to  be  thorough,  as  many  control  derivatives  as 

possible  were  implemented  in  the  aircraft  equations  of  motion.  Thus, 

Cp  was  derive'5  using  the  following  equation  and  the  data  is  in 
6h 


Table  III. 


4IH  .  i  .  ♦  ih  ‘X 
\  Six  hr  Xh 

h  AX 


(B-7) 


V/here  Cp  is  ;  .035  with  a  trim  angle  of  attack  at  *hj,  4.35°,  AX  is  the 

change  in  the  tail  incidence  angle,  and  Cp.  is  the  drag  of  the  unit 

*h 

horizontal  tail  (UHT)  at  the  incidence  angle  in  question.  The  average 

is  then  taken  on  the  summation  of  the  results  from  the  data  in  Table  III 

for  each  angle  change  to  determine  Cp  .  Using  the  average  applies  to 

°h 

the  derivation  of  all  the  remaining  control  derivatives  that  follow, 
except  for  C_  and  Cr,,.  . 


* 

v*** 


The  derivation  of  Cp  required  that  the  vertical  tail  be  modeled 

6r 
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ff 


as  a  wing  in  order  to  obtain  this  control  derivative  using  the  following 
formula . 


(B-8) 


Where  Cn  is  the  induced  drag  coefficient  and  Cn  ia  the  change 
ui  min 

in  minimum  drag  due  to  a  change  in  control  surface  deflection, 

Table  IV  contains  the  data  used  to  find  . 

4r 

For  Cq  change  in  drag  due  to  flap  deflection  was  found  using 
6f 

Equation  (B-8)  only  $r  is  now  The  data  used  in  the  derivation 
of  Cn,  ,  Ct ,  ,  and  (L,  is  in  Table  IV.  The  following  equations  were 


used  for  the  derivation  of  Ct  and  Cm  are: 

«f 


(B-9) 
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TABLE  III 


DATCOM  Data  for  Cn.  and  Cn. 

°h  ar 


°D«h 

\ 

5h 

CD 

«r 

^in 

% 

4.0°  Ref 

0.035 

mm 

0.00547 

9.57E-6 

6.0 

0.039 

BB 

0.00337 

3.88E-6 

8.0 

0.043 

-2.0 

0.00154 

6.917E-7 

10.0 

0.047 

0 

0 

0 

12.0 

0.050 

1 

0.00154 

1.8E-6 

15.0 

0.055 

0.00337 

6.1E-6 

1 

0.C0547 

1.29E-6 

*Units  deg"1 


TABLE  IV 


DATC  OM  Data  for  On.  ,  Cr  „  ,  and  C-, 

Of  Of  Of 


CD 

5r 

% 

sr 

C 

6f 

6f 

°Bmin 

(V 

Di 

CL 

C 

m 

5 

0.094 

-0.0137 

10 

0.00583 

6.19E-4 

0.187 

-0.0273 

15 

0.01356 

1.47E-3 

0.276 

-0.0402 

20 

0.02294 

2.59E-3 

0.361 

-0.0526 

25 

0.03453 

3.93E-3 

0.439 

-0.064 

30 

0.04743 

5.25E-3 

0.502 

-0.0732 

35 

0.06203 

6.34E-3 

0.549 

-0.0801 

40 

0.07578 

7.04E-3 

0.578 

-0.0843 

* Units  deg"1 


The  lateral  control  derivatives  ,  and  were  derived 

6h  $f 

with  the  data  in  Table  V  and  the  following  equations: 


(B-: 
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rS**NVH0ton*  »*fWl 


It-, 


The  control  derivative  was  calculated  using  the  following 
5h 

equation, 


n*  3 - 

°h  C.  .  2 

Xit  ( B-ll ) 


where  Ci  x  s  0.5937/rad,  Cn  =  0.2292/rad  from  (Ref  21)  end  Ci,  is 
xit  ui%  Ait 

twice  Ci  .  The  2  is  divide  through  to  obtain  the  magnitude  for  one 
5h 

control  surface  and  the  minus  sign  indicates  the  direction  of  yaw  with 
the  left  UHT  down. 


TABLE  V 


DATCOM  Data  for  Ci  ,  Cn,  ,  and  Ci , 
5h  6f  6f 


\ 

c 

n. 

5f 

% 

sr 

5h 

C1 

6f 

cn 

6.8 

3.329E-3 

5 

-1.228E-4 

2.8309E-3 

5.0 

2.466E-3 

10 

-2.456E-4 

5 . 6617E-3 

0 

0 

15 

-3.601E-4 

3.3004E-3 

7.5 

-3.6992-3 

20 

-4.193E-4 

9.6722E-3 

10.0 

-4.932E-3 

25 

-4.402E-4 

1.0157E-3 

12.5 

-6.164E-3 

30 

-4.791E-4 

1.1054E-3 

15.0 

-7.397E-3 

35 

-5.215E-4 

1.2031E-3 

17.5 

-8.630E-3 

40 

-5.646E-4 

1.3026E-3 

20.0 

-9.860E-3 

22.5 

-1.110E-2 

* Units  deg 
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TABLE  VI 
Data  for  C 


0. 

X 

Ss 

Cn  (DATCOM) 

Cn  (Ref  2) 

10 

3.47E-4 

3.0E-4 

20 

6.84E-4 

7.0E-4 

30 

1.0E-3 

1.8E-3 

40 

1.285E-3 

2.0E-3 

50 

1.53E-3 

2.5E-3 

60 

1.73E-3 

3.5E-3 

»Unit  deg 

To  obtain  a  degree  of  confidence  in  the  data  being  used,  a  corapar 
ison  was  made  between  the  A-7D  Aerodynamic  Data  (Ref  11:206)  and  the 
digital  DATCOM  data  for  CR  vs.  5^.  The  following  equation,  along  with 
the  data  from  Table  ”1  was  used  to  derive  Cn  . 


s  d. 


(B-12 


Since  the  data  of  (Ref  11 )  for  Gr  is  for  a  spoiler  deflector,  the 

average  valve  had  to  be  divided  by  a  scale  factor  of  1.75,  (Ref  14:62) 

in  order  to  obtain  a  Cnr  for  the  spoiler  alone.  The  values  are 

us 

1.91E- 3/rad  and  1.652E- 3/rad  for  the  A-7  Aerodynamic  and  DATCOM  data 


respectively,  with  a  difference  of  3-58E-4/rad.  This  is  an  acceptable 

difference  considering  the  different  methods  used  in  their  derivation. 

The  Cn  was  der.'red  with  the  following  equation 
6s 

^n.  •  b 
°DS 

3  OBS  (B-13) 

where  b  is  the  wing  span  and  OBS  is  the  out  board  end  of  the  spoiler 
which  is  43-465  of  the  semispan. 

This  completes  the  derivation  of  the  dimensional  control  derivative 

equations,  and  the  non-dimensional  control  derivatives  not  available  in 

(Ref  11)  or  (Ref  21).  After  reexamination  of  the  coefficient 

^r 

(which  may  be  zero),  it  is  suggested  that  this  control  derivative  be 
further  investigated  in  future  studies. 
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This  thesis  investigates  control  of  an  aircraft  when  there  is  a  primary  con¬ 
trol  surface  failure.  The  object  of  this  study  is  to  reconfigure  the  remaining 
control  surfaces  to  compensate  for  the  additional  forces  and  moments  generated 
by  the  inoperative  control  surface.  Tc  jtudy  this  flight  control  problem,  a 
comprehensive  aircraft  model  is  required  which  considers  each  control  surface 
operating  individually. 

A  six  degree-of-freedom  aircraft  model  is  developed,  including  all  the 
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individual  control  surfaces.  A  control  surface  input  can  produce  both  a  lateral 
and/or  a  longitudinal  response.  Thus,  the  equations  of  motion  cannot  be  de¬ 
coupled  for  the  design  of  the  control  laws.  The  coupling  between  the  axes 
requires  the  derivation  of  several  new  non-dimensional  control  derivatives. 

Using  the  geometrical  properties  of  the  aircraft  and  the  Digital  Datcom  computer 
program,  the  needed  control  derivatives  are  derived. 

With  a  comprehensive  aircraft  model  now  available,  the  entire  eigenstructure 
assignment  method  is  used  to  assign  both  the  eigenvalues  and  the  eigenvectors 
to  the  closed-loop  plant  matrix.  This  method  is  used  for  the  direct  digital 
design  of  a  multivariable  discrete  regulator  and  tracker  control  law.  The 
effect  of  increasing  the  number  of  control  inputs  on  the  relative  degree  of 
controllability  of  the  states  was  determined  by  singular  value  decomposition. 

This  thesis  concludes  that  a  direct  digital  design  for  reconfiguring  the 
multivariable  control  law  is  feasible.  However,  more  wind  tunnel  data  is 
essential  to  derive  the  additional  control  derivatives  for  a  more  accurate 
aircraft  model  driven  by  individual  control  surfaces.  Further  work  is  also 
necessary  to  perfect  the  assignment  of  the  closed-loop  eigenvalues  and 
eigenvectors. 
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